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Astigmatism of Skew Pencils in Optical Systems 
Containing Toric Surfaces’ 


Walter Weinstein? 


\ method is described for tracing skew rays through toric refracting surfaces and for 
determining the astigmatism of pencils around these rays. 


1. Introduction 


Toric refracting surfaces are used chiefly in two classes of optical systems, spectacle lenses, 
and anamorphotic photographic and projection objectives. In spectacles the toric surfaces 
are arranged to correct the axial astigmatism of the eye, and in anamorphotic systems two or 
more toric (possibly cylindrical) surfaces are included to give stigmatic imagery with different 
magnifications in different directions in the image plane. 

When the axial pencil is refracted by a toric,’ some astigmatism is introduced; if a pencil 
is incident obliquely, either centrally or eccentrically, astigmatism of different magnitude is 
introduced. It is therefore desirable in the case of both types of optical systems mentioned 
above to determine the astigmatism of obliquely incident pencils, in order to assess the imagery 
in the outer parts of the field. This is the problem which was studied in the work to be described. 

The problem arose in connection with spectacle lenses and may be outlined as follows: 

Figure 1 shows front and side views of a spectacle lens with toric surface; O is the center of 
the lens, C is the center of rotation of the eve, and P,, P2:, P; are points of incidence of the chief 
ravs of oblique pencils. Let the meridian and equatorial sections of the toric be, respectively, 
in the vertical and horizontal planes OP, and OP;,. Then the problem for chief rays passing 
through the axis and through ?; and 7 is trivial. The chief ray itself can be traced by the 
ordinary ray-tracing formulas for meridian rays and spherical surfaces, and the astigmatism 
can be determined by the usual methods for close sagittal and tangential fans, with the slight 
difference that the curvatures of the refracting surface in the two directions are not equal 


(see eq (3)). 





me 


FIGURE | 


If, however, the point of incidence of the chief ray does not lie in the meridian or equatorial 
section, as, for example, 7, in figure 1, the problem becomes more complicated. In the first 
place, even if the incident chief ray intersects the axis of the system, the refracted ray in general 
does not, so that the chief ray must for complete generality be considered to be skew; thus the 
first part of the problem is to trace a skew chief ray through the torie, 

lhe work described in this paper was Carried out while the author was a guest worker in the Optical Instruments Section, National Bureau 

standards, from January to July 1954 
Present address: Imperial College, London, England 
In what follows toric surface” will be abbreviated to “toric.” 
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Second, with oblique incidence as at any of the points /’;, 7’:, or Ps, astigmatism is produced 
in the refracted pencil for two reasons: (a) because the angle of incidence is not ZCTO, as in the 
ordinary case of spherical surfaces, and (b) because the refracting surface has two different 
principal curvatures at the point of incidence. In the case of Ps the principal directions of the 
astigmatism due to (a) are parallel and normal to the plane of incidence of the chief ray, i. e., 
approximately parallel and perpendicular to line, O7,; the principal directions of the astigma- 
tism due to (b) are, however, approximately parallel to OP, and OP,. Thus the resultant 
astigmatic pencil will have its focal lines in some intermediate directions, and the determina- 
tion of these directions further complicates the calculation It follows that the formal dis- 
tinction between sagittal and tangential directions can no longer be made; however, for the 
configurations of chief rey and toric that occur in spectacle lenses the rotation of the focal lines 
Is unlikely to be above 10°, so that there ts little diffieultv in relating the results for skew chief 
ravs to those for ravs in either principal section. 

The results of such calculations would constitute a survey of the astigmatism of a toric 
spectacle lens at all obliquities and azimuths. They would help to answer questions such as the 
following: What is the best bending of a toric spectacle lens? What range of powers can be 
manufactured with a common back curvature without departing appreciably from the best 
bending? Is there any advantage in making both surfaces of a spectacle lens toric in shape? 
Such questions as these arise in the design of so-called best-form spectacle lenses. These have 
previously been designed on the basis of computations in the two principal sections only 
However, an illustration is given by Henker ® showing the quality of an image produced by a 
pencil outside the principal sections of a typical sphero-toric lens of the Punktal best-form 
series; this illustration suggests that there is still room for improvement in the image. The 
work to be deseribed makes it possible to compute the imagery by skew pencils and so determine 
whether any improvements can in fact be made 

Section 2 is a recapitulation of the properties of astigmatic pencils that are required in 


deriving the astigmatism equations In section 3 the equations of a toric are given, and other 
properties are derived, in section | the ray -tracing equations for the skew chiel ray are derived 
and in section 5 the astigmatism formulas are derived Section 6 deals with the forms in which 
the input data must be given. In section 7 the mode of presentation of the results Is discussed, 


the results of a numerical example are given, and the possible uses of the computation are 


discussed 


2. Astigmatic Pencils 


Kor the purposes ol the present pauper i pene will he considered lo he made up ola narrow 
bundle of rays originating from a point in the object and surrounding a chief, or principal, ray, 
together with all the corresponding wave fronts The wave fronts are surfaces Of constant 
optical path length from the object pont and the ravs are normals to the wave fronts 

With a rectangular coordinate system with the « axis along the principal ray, the equation 


ol the wave tront passing through the origin will take the form 
r=Cyy? +eny ( terms of higher degree in y and 2; | 


linear terms are absent because the y plane must be tangent to the wave front at the origin 
Also it is possible to rotate the y and 2 axes about the + axis, so that the coefficient of the term 


in yz Vanishes when the equation takes the form 


I (CY c higher-degree terms. Z 
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Effects depending on terms in y and z of higher degree than the second will be neglected. 
To this degree of approximation it is easily shown that all rays of the pencil; i. e., normals to 
the surface (2) intersect one line through the point (l/e,, O, O) parallel to the 2 axis and another 
line through the point (I/e,, 0, O) parallel to the y axis, these being the two astigmatic focal 
lines. Thus this approximation takes account of astigmatism but not of aberrations such as 
spherical aberration and coma, which depend on higher powers of the aperture of the pencil. 
The coefficients have the significance that c, is the curvature of the section of the wave front by 
the sy plane and ¢,, that by the «2 plane 

Some results from differential geometry are now required. At any nonsingular point 7? 
of a surface the normal curvature (i. e., the curvature of the section of the surface by a plane 
through the normal) has a maximum and a minimum for two section planes which are at right 
angles to each other; furthermore, if these two planes, together with the tangent plane at 7’, are 
taken as coordinate planes, the equation of the surface in the neighborhood of 7? takes the form 
of eq (2), that is, the yz term is absent. The maximum and minimum normal curvatures are 
called the principal curvatures, and the corresponding section planes meet the tangent plane 
at /? in lines that lie along the principal directions at 77. Thus the principal directions are along 
the y and 2 axes for eq (2 
It is also known, by Euler’s theorem, that if a normal section plane makes an angle @ with 


the ry plane the curvature in this section Is 
Cla C. Gas & c sin a m 


Krom this it follows that these quantities are required to characterize an astigmatic 
pencil 
a) The direction of the chief ray ; this will be given by a unit vector, Q. with components 
A yr. ana Z 
b) The principal curvatures, ¢, and ¢,, of the wave front at a point on the chief ray 
¢) The principal directions on the wave front; these will be given by unit vectors R and S, 
and from the above it follows that QQ. R, and S are all perpendicular to each other. 


3. Toric Surface 


A toric surface is generated by the rotation of a circle (the generating circle, abbreviated 
to G. ©.) about an axis in its plane but not passing through its center, as in figure 2, a and 
b. In optical applications the lens surface is, of course, only a cap-shaped portion of the 
Whole surface, roughly centered on the pont A Such a surface is alined with a centered 
system of spherical surfaces by making the optical axis lie in the normal plane to the rotation 
axis, Which contains the center of the G. C. The optical axis then meets the toric in the 


pole \ 


ROTATION AXIS 


ROTATION AXIS 


G.C G.C 






OPTICAL A | OPTICAL [A | CENTER 
AXIS AXIS OF GC 


Q 
RING TORIC 


b 
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The toric refracting surface has different powers (in the optical sense) in different sections 
through the opticel axis: the powers are determined by the radius of the G. C. for the section 
in the plane of the rotation axis and by the distance from A to the rotation axis for that 
perpendicular to it; these are known as the meridian and equatorial sections, respectively. 

It can be seen that for given powers in the principal sections there are two possible forms 
of toric surface: the radius of the G. C. can be either less than or greater than the distance 
from A to the rotation axis, the two forms being known, respectively, as ring- and _ barrel- 
shaped torics. The ring form seems to be used almost exclusively for spectacles in commercial 
practice, but to the best of the author’s knowledge, this is merely because more rings than 
barrels of given curvatures can be surfaced in one operation, not for any optical advantage. 
The possibility of improvements from the use of the barrel shape seems not to be considered 


in the literature on “‘best-form” series. 
3.1. Equations of the Toric 


One can write down the equation of the toric surface by taking rectangular coordinate 
axes with the origin at the pole A of the toric refracting surface (fig. 3), the + axis along the 
optical axis, the y axis parallel to the rotation axis, and the ¢ axis in the equatorial plane. 
Let r, and r, be the radii of curvature in the ry and rz sections, respectively. Take angular 
parameters, 6 and @, as in the figure, so that @ is the longitude and @ the latitude of a point 


on the toric. Then the parametric equations of the toric are 


! } (7 r, +r, COs @) cos 6 


}) 


UV] r, SIN @ 
(7 r,+r, Cos @) sin i] 


It can be seen that eq (4) represent a ring- or a barrel-shaped toric according as 7, is less 
than or greater than r, in absolute magnitude,’ in both cases the rotation axis is parallel to the 
y axis. When the rotation axis is to be parallel to the 2 axis the equations are obtained be 
interchanging y and 2 throughout, including subscripts, in eq (4); this is true also for the 


equations in the remainder of section 3. 


ROTATION 
axis 


OPTICAL 
Axis 





FIGURE 3 


An explicit equation can be obtained by eliminating @ and @ in eq (4) and solving for 4 


r=Tf, 7 / r,*—y*)'2}-— 2° +" 5) 


The ambiguous signs are due to the fact that there are four (possibly complex) values of 
x for any given y and 2, since the toric is a quartic surface; the required solution is chat which 
tends to zero as y and 2 tend to zero, so that the correct signs are, respectively, minus and 
plus. When these are inserted the result is unsuitable for numerical computation when y and 





’ Provided r, and r, have the same sign, which is generally the case in spectacle lenses 
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zare smell compared tor, and 7,, as is generally the case. This difficulty can be circumvented 
by transforming the equation by means of conjugate surds in the usual way and using curva- 
tures instead of radii, giving 


G) G- |G" c.F]! 
where 
F=¢,2@¢ + 2¢,7G—c,7%c.y', G=1+ (1 c,*y")'?. 


In this form the equation of the toric is suitable for numerical computation with all 
relevant ranges of the variables. 


3.2. Principal Directions and Curvatures 


In considering the refraction of an astigmatic pencil through a toric, it is necessary to know 
the shape of the toric in the neighborhood of the point of incidence. Thus, as in section 2, the 
following are required: 

(a) The direction of the normal at the point of incidence; this will be given by a unit 
vector, U, with components, A, uw, and »v. 

(b) The principal curvatures, c, and ¢,, at the point of incidence. 

(c) The principal directions at the point of incidence; these will be specified by unit vectors, 
V and W, and again the vectors, U, V, and W, form a trirectangular system. These quantities 
are easily determined for a current point (z,y,2); if the toric has principal curvatures, ¢c, and ¢,, 
at the origin and if the rotation axis is parallel to the y axis, as in section 2.1, the results are 


[° (cOs @ COS 6, sin 2, Cos @ sIn ), 7) 
, ° Cc. . 
o=—aresin ¢,y,¢6—aresin —_, ? 
; G—c.lyy 
c c 
Ge COs @O \) 
Cc . 
ee 
4 sin o COS 4, COS o, sin o sin @ 
y 
I} = (sin @, 0), cos @ J 


4. Tracing the Skew Chief Ray 


Using the quantities defined in sections 2 and 3.2, the problem can be stated as follows: 

Given @ and the coordinates of a point on the chief ray, the coordinates (z,y,2) of the 
point of incidence of the ray on the toric must be determined (transfer process); next U atthe 
point of incidence must be found and Q’, the direction of the refracted chief ray, determined 
(refraction process). The quantities c,, ¢,, R, and S will be known for the incident pencil 
and ¢,, x, V, and W can be found, since the point of incidence is known; from these the 
quantities ¢), ¢;, R’, and S’ for the refracted pencil must be determined. This last, the astig- 
matism calculation, will be described in section 5. 


4.1. Transfer Process 


The transfer problem is that of finding the intersection of a straight line with a toric. A 
direct approach would require the solution of an algebraic equation of the fourth degree; this 
would be a laborious process, so an iterative method was adopted. In addition to the general 
alvantage of iterative methods, that mistakes in the early stages of the calculation are auto- 
matically eliminated at the cost of extra computing time, this method also selects the correct 
solution from the four possible solutions. 
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Let P,; be a point on the ray near the required point of incidence.S A point 7?) on the 
toric is found with the same y and 2 coordinates and the tangent plane 7) at 7’; is determined 
The next approximation 77, Is taken at the intersection of the ray with 7) and the 
evele P Me -T.— PRP; . repeated until the coordinates of ?, and ?’,,., are the same to the 
required number of significant figures. These are then the coordinates of 7’, the point of 


ine iden ce 


Let the coordinates of P, be (4, .¥,, and those of P. be (2).y,.2,); then ri is obtained 
by substituting y, and in eq (6 The equation of the plane 7), is A, (4—v., u, (YY, 
0), 
where \,, u,, and vz are given by eq (7) so that 4 is given by 
/ / ) 


AN (sr, I 10 


\lso 
y. 
/ Y,, 4 ¥° 
Z 1 
* 
This completes the iterative evele It can be shown to be a second-ordet iteration,” that Is, 


is of order of magnitude 62, so that it is always rapidly convergent, provided the initial 
approximation is good enough 
The choice of the initial approximation depends on the available means of computation 


If a high-speed automatic-computing engine is available, 7; mav be taken as the point where 


the ray meets the yz plane. Let (a %, be the coordinates of the point 7?) where the ray 
leaves the previous refracting surface (or the entrance pupil then 7, is given by 
0) 
Y] / rod A Z 
AA 


Starting from this approximation, some four or five iterations would probably be meeded to 
reach a result correct to six significant figures 

If a desk computing machine is to be used, an appreciable saving in computing time can be 
effected by using an improved first estimate. It usually happens that chief rays directed at 
the same pupil have already been traced in the meridian and equatorial sections of the toric 
Let the « coordinates of their points of incidence be «, and s,, respectively The section of 
the toric by a plane through the # axis and making an angle @ with the ry plane Is approximately 
a circle, of curvature ¢,cos*? a+c¢, sin’ a, so that a chief ray in this plane will meet the torie at a 


point of which the + coordinate is approximately «, cos* a+r, sin’? a This P; is given by 


J r. COS” yr, sill a 


This approximation is usually so good that only one iteration is needed to obtain six 


, 


significant figures in the coordinates of / 
Weinstein, I Phy : London) (B) 63 
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4.2. Retraction Process 


The refraction calculation for the chief ray can now be carried out by well-known formulas 
Let J be the angle of incidence and n the index of refraction of the medium in which the ray is 
incident; let primes denote corresponding refracted quantities. Then 


cos 1=@Q.-U, 14 
cos [’ t ( | = 1—cos* / ) * (15) 
Q’—~, Q+( cos 1’—~, cos 1) U Lt 


In eq (14) Q-U denotes the scalar product of Q@ and U 
This completes the refraction process and the coordinates ana direction cosines of the 
refracted ray are then used, with a suitable change of origin, for the incident ray at the next 


surface 


5. Astigmatism Calculation 


In order to determine equations giving R’, S’, ¢{, and ¢), in terms of Known quantities, a 
eeneral theorem is required Referring to figure 4, let 7? be the pomt of incidence of the chief 
ray, and let 2 be the pomt of incidence of another ray of the pencil Suppose the ineident and 
refracted wavefronts through 7? to be constructed, and let the ray through B meet these in ¢ 


ana ES; respectively Then since the optical path length between two wavetronts ts the same 








along any ray, it follows that n»OB—n’C"B’ or, if CB is denoted by 7 
A nq) (). 


Here the svmbol, A, means that the refracted quantity minus the meident CpUlaari tits is to be 
taken. This relationship is true for any relative positions of 7? and 6B and for any shape of 
refracting surtace 

In order to use eg 17) for the present purpose, ¢ ania q will be expressed in terms of the 
coordinates of the point Boon the refracting surface and the parameters of the wavefronts 
The resulting equation will give the normal curvature of the refracted wavefront in a section 
containing the ray through Band by differentiating with respect to the azimuth of 6 the maxi 
mum and minimum normal curvatures, ¢ and ¢!, can be found. The same calculation gives 


the directions R’ and S’ of the principal SCCTIONS, 


INCIOENT REFRACTED REFRACTI 
NG 
WAVEFRONT 
WAVEFRONT SURFACE 


- 

(n) 
Q - 4 
CHIEF RAY 
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5.1. Normal Curvature of the Refracted Wavefront 


Figure 5 shows the refracting surface near the point of incidence. Let 7? and B have the 
same significance as in figure 4 and let Vand W be the unit vectors along the principal directions. 
Let B, be the projection of B on the VW plane; let a be the angle between PB, and V, and 
let PB be denoted by p. Then the position of B relative to P will be given by p and a@ as 


coordinates 





In order to express g in terms of p and a, use is made of the fact that CB in figure 4 is a 


normal to the wavefront It is found see section 5 la that 


q Pp Q.V cosa Q-W sina +p qQ Ule COs” a-~-Cc, SID a c R V cosa ’ R W sina 


c,(8-V cos a+ S-W sin a)* terms in Pp. IS 
and there is a similar equation for q’. 
Substituting for g and q’ in eq (17 
pAn Q-V cosa ! Q-W sin ain Sp An Q Uie, cos* as Cy sin’ a c(R-V cos a R.W sine 
¢(S-V cos a+S-W sin a)*} +-terms in Pp 0). 1%) 


Only terms up to and including the square of p are taken since, as explained in section 2, 
this degree of approximation takes account of astigmatism. 

Equation 19 holds good for any values of p and a, so that the coefficients of different 
powers of p may be separately equated to zero. If this is done, the linear term vields an altern- 


ative form of Snell’s law and on putting a=0 and 7/2, the equations 


nQ-W|—0 20) 





Alng@-V } A 


are obtained. These could serve as a check on the refraction calculation, 
Equating the coefficient of p? in eq (19) to zero and dividing through by 1—(Q’ V cos a 


Q’-W sin «)* gives 


( R’.V COs an R’.W Sin a) Cc. Ss’ .7 COS a S’.W sin a) 
| Q’.V cos a+Q’-W sin a)’ 


iin R.V COS an R.W Sill ald He S.V COS a S- Ww Sih a) CC. cos” a Cc sin a An@ U 2] 
n’ (1—(Q’-V cos a+ Q’-W sin a)’ 7 


Now the equation of the refracted wavefront can be shown to be 


c(R’.V cos at R’-W sin a)?+c)(S’-V cos at+S’-W sin a)*|(a° + y 
R’.V cos a~ R’.W sin ay- S’.V Cos an S’oW sin a) 


taking the z, y, and ¢ axes along Q’, R’, and 8S’. Thus the left hand side of eq (21 is ec’ (a), the 
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normal curvature of the refracted wavefront in a section containing the point B with angular 


coordinate @ on the toric. 
a. Derivation of Equation (18) 


Two auxiliary coordinate systems are required, each with origin at 7, the point of inci- 
dence. One, the x, y, 2 system, has its axes respectively parallel to Q, R, and S; the other, 
‘, y’, and 2’ has its axes parallel to U, V, and W. The two are related by the equations 


r= (Q-U)x’ + (Q-V)y’+(Q-W)2’, y=(R-U)r’ +(R-V)y’'+(R-W 
>= (§-U)x’ +(S-V)y’ +(S-W 


If the point Bon the refracting surface (fie. 4) has polar coordinates p, a), its coordinates 


in the 2’, y’, 2 svstem are 
“ h(c, cos* a+¢, sim a P+O(py), youpcosarO(p), ‘<=p sin a+O(p 
Let B have coordinates (uw, 7, w) in the (x, y, svstem; then 


u—3(Q-U)(c, cos’ a+e, sin® a)p Q.V)p cos a Q-W)p sin a+ Op 
R-V)p cos a+ (R-W)p sin at O(p 
w— (S-V)p cos a+ (S-W)p sin a+ O(p 
Let (4, y be the coordinates of Cin the wv, y, 2 svstem, so that 
r=3(c¢,y’+c Oy, 


The direction cosines of the normal to the wavefront at C are, to the approximation O(y 


& CM. ( and since ('B is a rav, or wavefront normal, it follows that 
ul J q, Y CYY, rn Co 
There are now enough relations to eliminate vs, y, 2, vu, ¢, and w. The resulting equation involy- 


Inge p, g, ana a can be solved for , remembering that terms W hich are O Pd cunh be neglected, 


since only astigmatism is in question. The solution for g is found to be eq (18). 


5.2. The Principal Curvatures of the Refracted Wavefront 


The principal curvatures ¢) and ¢) are the maximum and minimum values of e¢’(a 
Since all quantities on the right of eq (21) are known with the exception of a, ¢) and ¢! can be 


obtained in the usual with by differentiating with respect to a ana equaling to zero Let 


a —3{ne,((R-V)?-+ (R-W)*) + ne, (S-V)?-+ (SW ¢, +¢,)dnQ-U], 2» 
b—3[ne,((R-V)?— (R-W)?) + ne,((S-V)?—(S-W c,+e,)AnQ@-U], 23 
c—ne,(R-V)(R-W) +ne,(S- V)(S-W), 24 
d—3((2—(Q’-V)?—(Q’-W)?), 25 
e—1((Q’.V)?—(Q’-W)?), oF 
Q.ViiQ ew 27 


The equation for @ ts 
ea Ad sIn- a de at) Cos” a hf ai 








and the two solutions are 


y 3 r Y re e 
7 » a > we -S 
where 
se at . 
oO—aresin ov 
ea Ad de at 
ana 
Af é¢ 
Y aresin . 30 
ea Ad de at 


The values of ¢/ and ¢! are then obtained by putting a 


a, and a in the right hand side of 
eq (21 


5.3. Principal Directions on the Refracted Wavefront 


In figure 6, let 7? be the point of incidence, ?’B the refracting surface, and P¢ 
wavefront Let Q’ be the unit vector along the chief ray and let (’°B bi 
astigmatic pencil Let the projections of ¢ 


‘che refracted 
‘another ray of the 
‘and / on the tangent planes to the wavefront and 
refracting surface be (, and B, respectively. Then V cos a 


W sin ais a vector in the direction 
PB _so that a unit vector in the direction re 


‘is given by 


V cos a W sin a qQ’ qQ’ V cos a-t QW Sill a 7 
1—(Q’-V cos a+ Q’-W sin a)*|"” . : 


Thus, the vectors R’ and S’ are obtained by putting a=a, and a in the expression (31 It 


should be noticed that in this derivation it was assumed that the points ¢ 


‘and B hie On a plane 
through Q’, 1. e., that the rav through CC’ intersects the chief ray. 


This is not strictly true in 
veneral, but the errol depends on terms 1n the equation of the wavefront of the third degree 1T) 
the aperture, so that astigmatism calculations are not affected 


REFRACTED Vcosa+Wsina 

WAVEFRONT 

REFRACTING 
SURFACE 





Ts) 


This completes the retraction computation for astigmatism The zx and S’ vectors are 


transferred without change to be the new R and S vectors at the next surface 
transter the curvatures the distance L 


is easily found to be (#—+,) VY. in the notation of section 4.1 Then the new curvatures are 


In order to 
between surfaces along the chief ray is required; this 


given by 


ls l/e,)’—L ) - 
1/e,.)’—1 , 7 
( { sJ 
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6. Input Data 


If the wavefront incident on the first surface is not already astigmatic, some arbitrary 
directions for R and S must be assigned. Similarly if a refracting surface is spherical, dire 
tions for Vand W must be assigned 

These pairs of vectors must form trirectangular sets with @ and U, respectively, but subject 
to this condition they can be chosen to suit computational convenience. Choosing R to be 


coplanar with @ and the y axis, it is found that 


R AY 1-—¥?*,-—Y¥Z 


2,0,—A 


S : 
l—} 


Since )’: Is rarely vreater than 0.1, the Vectors so defined ure well determined numerically 


In a similar wav V and W ean be taken as 


, Au, | Me, LV 
J 
is 
0,—X - 
, | / 
WW 
l—yp 
It is also necessary to choose between rine and barre] shapes for the Lol ic surfaces, In the 
equations of section 3.1 the rotation axis of the toric is taken parallel to the y axis and the toric 
will be ring or barrel shape according as c, is numerically geater or less than ¢,. Thus, if it is 
required to have a barrel shape with ¢, c, or a ring shape with (¢, c,, the equations 


must be changed so that the rotation axis is parallel to the 2 axis. This can be done by inter- 
changing y and in the symbols and subscripts of eq (4) to (8), inclusive, as mentioned in 
section 3.1 


7. Presentation of Results 


7.1. A Numerical Example 


In the application to spectacle lenses the results are presented in the form of a diagram 
giving the astigmatism of a number of pencils and the directions of the principal sections on 


the wave front. An example is given in figure 7, 


5933 
3 92! 





RING TORIC BARREL TORIC 
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The wavefront curvatures are marked on the appropriate sections, in diopters with respect 
to the vertex sphere. The apparent rotation of the principal sections of the wavefront can be 
found from the consideration that the vertical plane containing the chief ray contains the 
direction which is apparently vertical to the eve. This plane is the plane through @Q and 


parallel to the y axis. From this it can be shown that the angle, A, through which the R and 


S vectors appear to be turned is given by 
tan 1 ’ (3) 


where ry and s, are the y components of R and S. 
The nominal powers of the lens in this example were 5.95 diopters in the vertical section 
and 3.62 diopters in the horizontal section The curvatures of the front (toric) surface were 
0.2199 em”! and + 0.1839 em™', that of the back (spherical) surface was —0.1147 em~', the 
center thickness was 0.48 em and the refractive index was taken as 1.523. The chief ray was 
made to intersect the axis 27 mm behind the lens at an angle of 30°. The results for both ring 
and barrel shapes with the same nominal curvatures are given. 
Computations have also been done on the prescriptions of a number of lenses and measure- 
ments have been made on these with a vertex power measuring instrument. The agreement 


is satisfactory within the limits of experimental error. 


7.2. Discussion 


In the numerical example given above the rotation of the focal lines in intermediate 
azimuths is negligible; this may be expected to be the case for all but very heavily cvlindrical 
spectacle lenses of best-form type. It can also be seen that the variation of astigmatism over 
the field is less for the barrel than for the ring shape (0.14 D and 0.21 D, respectively). This, 
of course, is not a general conclusion, but it suggests that further computation would be 
desirable. It is possible that for heavy evlinders a considerable improvement might be 
obtained by using barrel shapes, but this might be outweighed by the increased manufacturing 
costs. This would depend also on visual tolerances for astigmatism; it seems that more 
experimental work is required on this subject 

A different question is that of the wecuracy W ith which tories are produced The principal 
curvatures check within the limits of accuracy of a lens gage, and the sagittal curvatures change 
or remain constant along the meridan or equatorial sections, respectively, to the same accuracy 
However, this may not be enough to prove that the lens surface is a true toric to an accuracy 
which justifies the calculations for chief ravs in intermediate azimuths. The mode of manu- 
facture is certainly such as to produce a very close approximation to a torte if the surfacing 
machine really acts in the intended manner, but it is quite possible that distortions and non- 


No 


measurements other than checks on the principal curvatures seem to have been made A 


uniform wear in the tools and machine parts may cause departures from the toric shape 
detailed check on the « coordinates of a network of pomts would be needed 

Thanks are due to my friends and colleagues of the Optical Instruments Section for thei 
interest in this work, in particular O. N. Stavroudis and Lovd Sutton; the latter did the 
numerical computations and checked the mathematics 


WASHINGTON, June 29, 1954. 
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Osculatory Interpolation in the Complex Plane 
Herbert E. Salzer' 


Tables of coefficients to facilitate osculatory n-point interpolation (m= 2(1)7) in the complex 


plane are given. 


The writer has described in a previous article 





1} 2a new method for osculatory interpolation 
that can be applied also in the complex plane when one has an analytic function f(2), 2=2+-7y, 
which is known for arguments 2 at equal intervals along any straight line in the 2-plane 
However, for the tabular arguments 2,=.,+/7y, in the form of a Cartesian grid of length A, 
greater accuracy may be had by basing the Hermite osculatory interpolation formula upon 
values of f(2,) and f’(2,) at points =, which are closer together and not necessarily lying upon 
a straight line (except, of course, in the 2-point case). For a detailed discussion of the Hermite 
osculatory interpolation formula, the reader is referred to the previous article [1]. This 
present paper merely aims to supplement that article which was intended primarily for real 
functions, by giving here the corresponding auxiliary quantities, a, and &, which are suited 
better for complex interpolation. Whereas the quantities, a, and ), for real interpolation 
were tabulated up to the 11-point case (21st-degree accuracy), this present tabulation for 
complex interpolation does not go beyond the 7-point case (13th-degree accuracy), which is 
more than adequate for most of the practical problems that would arise. 

To interpolate for f(z), where z=2)+Ph and P=p+iq is now complex, we choose the 

~y +kh for k equal to certain small complex integers and such that 2) will always be the 
lower left corner point of the configuration of points, 2) and 2,. Also, in almost every example, 
p <1, ¢ <1, so that the point = is usually not outside the square whose corners are 2, 2), 24 
and 2, For each of the n-point formulas, n=2(1)7, the points 2, are chosen to lie im the 
following configurations: 


Two-point Three-point four-point Five-point 
t ~1 ! i 
1 
Six-pomt Seven-point 
2144 
1 1 c . 
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Hermite’s n-point osculatory interpolation formula for /(z)=f, in terms of f(- f. and 
/’,, is expressible in the concise form 


where 
x, =a, (P—k 6,/(P—| 2 

wna 
3, =a,./(P—k ; 
the summation being taken for /£ corresponding to the n-points of the configuration. The 


dependence of ay, By, a, and 6, upon nm is understood. The auxiliary quantities, a, and by, 
are tabulated below for »=2(1)7 for each of the above n-point configurations of points 


The actual definitions of a. and 4, are as follows: 


at 


h 2C(n ppl / J } | {1 } ) . (5) 


where the product is taken for ) corresponding to the n—1 points 7,72, of the configuration 
and where for each n, C(n) is chosen as the least integer which allows the quantities, a, and 4,, 


tO appear as complex integers 


lal ! } 
wo-point SiX-point 
| / 2 a 500 h 2000 — 2000 
| / 2 a 1280 — 960 I, IS56 GB2BOS 
a 30 10 hy 83 + 269 
a 1280 — 960 h, H20O8 1856 
Three-point a 2000 /, HO000 ~— 6000 
a 30— 40 } 2690 + 83 
2 } } 
| 3 
3 Seve point 
SO 60 h 648 + 64 
Four-po nt a ISO H40 hy 3264 P0418 
30 10) 2 3 324 
l 2 } 320 — 240 h 1872 - 704 
; ; a 2000 by, HO0OO — 2000 
} }, 3 2 st) ht 240 14S 
i, 3 3 a 1O ho, 27 57 
Five 
125 } 375 500 
a S00 } 1500 500 
20+ 15 117 14 
SO. 60 } 50S +4 
250 } S00 O00 
Formulas (1), (2), and (3) for complex osculatory interpolation, which utilize the Cartesian 


eri, will be espe ially useful In connection with such tables for complex arguments as | 


logarithm of the gamma function. log I] together with its derivative the psi function, ¥ 
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Je , 4 and 


2) Bessel functions of the first or second kind giving J and ./, 


e-“ du with its 


(2 Ke _ or linear combinations of them, (3 probability Integra 


integrand, (4) miscellaneous tables of integrals of the more clementary functions where the 


first derivative or the integrand, although not tabulated, is easy to calculate, namely, the fune- 
Hon e~“/u)du and (5) tables of solutions of important linear differential equations, together 
* 

with their first derivatives [2]. In all such tables, and in many others, the user will find these 
complex osculatory interpolation formulas to be particularly convenient. ‘They are especially 
suitable in those cases where the grid length / is too large for sufficiently accurate complex 


(3, 4) or formulas cor 


interpolation, using either tables of complex interpolation coefficients 


)) above [5], where the tables or formulas involve the use of only 


responding to (1), (2), and (3) 
the functional values / 
Illustration.—_To demonstrate the use of these tables in formulas (1) 
an example from [2], where the modified Hankel functions Aj(2)=(4 2)77.°(4 2), ped, 2, 
0.1 in the complex plane. 


to (4), consider 


and 


their first derivatives are tabulated over a Cartesian grid of length / 


Suppose that it is required to find A, (1.24579 316+4-0.96155 8037), using the four-point 
osculatory interpolation formula and the tabulated values of hy, Ay at 2y9= 1.24-0.97, 2 = 1.34-0.90, 
1.2+ 1.02, 1341.07 [2, pp. 21, 23]. Employing (2), (3) and a,, 6; from these tables, 


since P—0.45793 16-+-0.61558 O37. one finds that 


a (0.45793 164-0.61558 037)" 3-+-37)/(0.45793 16+-0.61558 O22), 
a i 0.54206 84+-0.61558 037)?+ (—38+32)/(—0.54206 844-0.61558 031), 
ae //(0.45793 16—0.38441 977)?+ (8—37)/(0.45793 16—0.38441 977), 
a 0.54206 84—0.38441 977)" 3—37)/(—0.54206 84—0.38441 977), 
B, = + (0.45793 164-0.61558 037), B; i/(—0.54206 84+-0.61558 037), 
B (/(0.45793 16—0.38441 977), B i/(—0.54206 84—O0.38441 977), 


from which one obtains 


ay = 7.09824 5538— 1.29187 90237, By = 1.04576 6212+-0.77794 TSS89:, 
a 6.63654 O0823--0.51581 36957. B, 0.91498 3814-+-0.80571 74872, 
a, ~ 9.82415 3465—1.10147 8051), B= 1.07536 0453 — 1.28099 97331, 
a 8.43081 1380+-1.81986 93857, re} 0.87048 2793 — 1.22746 386677, 
Then for (1) we employ the following tabulated values of A, and h; 
f,— 0.19018 585—0.19313 8407, fy ~ 0.14952 6144-0.35229 4614, 
f—0.20311 754—0.15966 4187. f,—O.10901 235+-0.33659 2310, 
f,—0.15891 O89—0.17847 850i, f, 0.14376 041 -+-0.29875 2520, 
f 0.17143 575—0.14888 10S; fy 0.10669 347 --0.29764 346, 
and also 
Say, = 31.98975 121 —0.05767 399) 
so that we have 
hy(1.24579 3164-0.96155 803%) =f~lahi+ af, +aif 
ji a oil 0.1! Bofi+BihitBilitBishies }V/SSe 
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whose terms are 


[(7.09824 5538— 1.29187 90237)(0.19018 585—0.193813 8407 

(6.63654 0823 +-0.51581 36957)(0.20311 754—0.15966 4187 
(9.82415 3465—1.10147 80517)(0.15891 O89—0.17847 850i 
8.43081 1380+ 1.81986 93857)(0.17143 575—0.14888 1087) 

-O.1. (1.04576 6212+-0.77794 78897)(0.14952 6144-0.33229 4617 
(—0.91498 3814+0.80571 74877)(0.10901 235-4-0.33659 2317) 
(1.07536 0453—1.28099 97337)(0.14376 041+0.29375 2327) 
(—O.87048 2793 — 1.22746 36677)(0.10669 347 +-0.29764 3467) >| 


31.98975 121—0.05767 3997). 


or after multiplying, 


1.10047 4414— 1.61664 O8967) + (1.43035 4817—0.95484 70407 
(1.36457 4820—1.92843 70327) + (1.71628 6592—0.94319 76317) 
0.1 0.10213 8505-+-0.46382 60217) + (—0.37094 2846—0.22014 33597 
0.53089 0903 +0.13173 25817) + (0.27247 1703—0.39005 58687) °] 


31.98975 121—0.05767 399) 
5.64471 8767—5.44458 66627)/(31.98975 121—0.05767 3997) 


or finally, 

0.17676 025—0.16987 9167, 
which is correct to within a unit in the last place, as was seen from the independent calculation 
of A,(1.24579 316+-0.96155 8037) from the Taylor series arouad 2 — 1.2 +-0.97. 


1] H. KE. Salzer, New formulas for facilitating osculatory interpolation, J. Research NBS §2, 211 (1954 
RP2491 

2} Harvard Computation Laboratory, Tables of the modified Hankel functions of order one-third and of 
their derivatives (Harvard University Press, Cambridge, Mass., 1945 

3] A. N. Lowan and H. E. Salzer, Coefficients for interpolation within a square grid in the complex plane 
J. Math. Phys. 23, 156 (1944 

1} H. E. Salzer, Coefficients for complex quartic, quintic, and sextie interpolation within a square grid, J 
Math. Phys. 27, 136 (1948 


>} H. FE. Salzer, Formulas for complex Cartesian interpolation of higher degree, J. Math. Phys. 28, 200 (1949 
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Macromethods For Reducing Aldonic Lactones To Sugars: 


The Preparation of L-Glucose’ 


Harriet L. Frush and Horace S. Isbell 


Two macroprocedures are described for the efficient reduction of lactones of aldonic 


acids by sodium amalgam in the presence of an oxalate buffer. 
flask equipped with a stirrer that sweeps the bottom of the flask 


‘“‘blender.”’ 


quantities. By the procedures described, 


One employs a stainless-steel 
; the other, a high-speed 


Either method is suitable for the preparation of rare sugars in relatively large 
L-glucose 


has been prepared from .L-glucono-é- 


lactone in 88-percent yield, for use as carrier in the preparation of L-glucose-1-C', 


1. Discussion 


Synthesis of sugars by the evanohydrin method 
of Kiliani and Fischer has been, and still remains, 
of great importance to carbohydrate chemistry. 
Originally, the process was tedious and the vields 
were low. Thus, in the classical synthesis of L- 
vlucose from L-arabinose the over-all vield was only 

percent [1].2. Although the technique of reduction 
by sodium amalgam has since been improved [2, 
1], the procedure has nevertheless been laborious 
and somewhat variable as to results. For use in the 
synthesis of C'*-labeled sugars on a semimicroscale, 
a convenient procedure was recently developed for 
the sodium-amalgam reduction of millimole quan- 
tities of aldonic lactones in vields as high as 85 to 
90 percent [5, 6]. The reduction is carried out in a 
closed system. Accurate pH control is obtained by 
using an excess of a slightly soluble acid salt (sodium 
acid oxalate) as a buffering agent. This avoids the 
necessity for continual adjustment of pH by titration 
of the alkali formed in the reaction. The method 
has now been adapted to macropreparations with 


equally successful results. Two procedures were 
found to be convenient for the reduction of a 5- or 
10-¢ quantity of lactone. The simpler and more 


rapid procedure involves treatment of the lactone 


with sodium amalgam and = excess sodium acid 
oxalate in a high-speed ‘blender’ partially filled 
with ice. The maximum vield of sugar from the 


lactone compares favor ably with that obtained in 
the semimicro apparatus [7], and requires approxi- 
mately the same proportion of amalgam; the reaction 
is complete in 15 min. For the preparation of C'*- 
labeled sugars, however, the procedure is complicated 
by the formation of mist, or spray, that is difficult 
to confine to the apparatus. The second procedure 
for reduction makes use of a stainless-steel flask 
equipped with a stirrer that sweeps the bottom of 
the flask. The amount of sodium amalgam necessary 
for maximum vield of the sugar was found to depend 
somewhat on the efficiency of stirring, and, in the 


apparatus emploved, was 2 or 3 times as large per 





thesis 


mole of lactone as that required in the blender. 
However, the two procedures give comparable 
results, and both are suitable for the macroreduction 
of any aldonic lactone.* 


2. Experimental Procedures 


2.1. Reduction of L-Glucono-é-Lactone to L-Glucose 
in a High-Speed Blender 

Approximately 500 ml of ice slush was agitated 
with 55 of sodium acid oxalate and 10 ¢ of 
L-glucono-d-lactone in a high-speed blender having 
stainless-steel blades. After a few seconds, 130 g of 
5-percent beter yen amalgam pellets [8] was added 
quickly; agitation was continued for 15 min, during 
which time the temperature rose to 33° C, The 
mercury was then separated, and the solution was 
neutralized with dilute sodium hvdroxide solution 
until a faint but permanent pink color was obtained 
in the presence of phenolphthalein indicator. The 
solution was evaporated under reduced pressure to a 
volume of about 100 ml and treated with 5 volumes 
of methanol. The precipitated salts were separated, 
washed with a little methanol, and discarded. The 
filtrate was concentrated under reduced pressure to 
about 50 ml, and again treated with 5 volumes of 
methanol, Precipitated salts were removed by fil- 
tration, and the solution, after concentration under 
reduced pressure to about 50 ml, was deionized by 
passage through 60 ml of mixed cation and anion 
exchange resins in a column having a laver of cation 
exchange resin at the bottom. The neutral solution 
and washings from the column were evaporated un- 
der reduced pressure to a thin sirup (about 10 ml) 
which was diluted with an equal amount of methanol, 
and then with 2-propanol almost to the point of 
turbidity. The solution gave 8.92 g of crystalline 
L-glucose in three crops, corresponding to 88.3 percent 
of the theoretical vield. 


or 

= 
e 
0- 


The aldonic lactones vary considerably in ease of reduction and in the tend 
v of the sugar to be reduced further to the corresponding sugar alcohol Con 
quently e optimum amount of amalgam to be used must be determined for 

h lactor iu i 
Equival quantities of Amberlite 1R-120-H, analytical grade, Rohm A 
Haas Co., Philadelphia, Pa., and Duolite A-4, Chemical Process Co., Redwood 
Cit Cal The laver of cation exchange resin at the bottom of the columr 
i he effluent when insufficient mixed resir used 








2.2. Reduction of D-Glucono-é-Lactone-1-C'' on a 
Relatively Large Scale in a Stainless-Steel Flask 


Mighty grams of sodium acid oxalate, 7 g of crude 
p-glucono-é-lactone-1-C™ having an activity of 5,330 
uc, and 500 ml of ice water were placed in a two- 
necked, 1-liter, stainless-steel flask equipped with an 
efficient stirrer” The solution was stirred vigorously 
and ISO g of 5-percent sodium amalgam pellets was 
added at one time. After 2 hr of stirring at ice-bath 
temperature, the mercury was separated ; the reaction 
mixture was made alkaline to phenolphthalein and 
then carefully neutralized with oxalic acid.° The 
crystalline sodium salts were separated and washed 
with methanol. The combined aqueous solution and 
methanol wash liquor were concentrated under re- 
duced pressure to about 250 ml, after which a crop 
of crystalline sodium oxalate was separated by filtra- 
tion, washed with methanol, and discarded. The 
filtrate was concentrated to about 20 ml at a temper- 
ature less than 40° C, and another crop of salts 
was separated. The resulting solution was concen- 
trated under reduced pressure at 30° C to a sirup 
that was then diluted with 5 volumes of methanol 
The mixture was filtered, and the residue was washed 
with 30 ml of methanol in 3 portions. Finally, the 
methanol solution was diluted with an equal volume 
of ethanol, and the salts that precipitated were sep- 
arated by filtration and washed with methanol. The 
filtrate was concentrated to remove most of the al- 
cohol; the resulting solution was diluted with 50 ml 
of water and was passed over a column containing 100 


ml of mixed cation and anion exchange resins. The 
| t T i simle steel rod fitted with a 
r I the fla k vd Vel t } 
| } kept neutral during t | ‘ ' 
! itl 


column was washed with water, the combined solu- 
tion‘ and washings were evaporated to 50 ml under 
reduced pressure at a temperature less than 40° C 
and the solution was then freeze-dried. The residue 
Was taken up at oncein 15 mlof methanol and brought 
to crystallization by the addition of 2-propanol to the 
point of incipient turbidity. The crystalline p-glu- 
cose-1-C that formed, after separation and reerystal- 
lization, weighed 5.47 ¢ and contained 3.954 uc of 
carbon-14. By use of nonradioactive carrier, an ad- 
ditional 658 we of crystalline a-p-glucose-1-C'’ was 
obtained. Thus the radiochemical vield was 4,612 
uc, or 84.3 percent of the theoretical.” By elution 
with 10-percent aqueous acetic acid of the mixed 
resin used for the purification of the sugar, crude p- 
gluconic acid-1-C" with an activity of 645 we was 
recovered, 


” The Hl of the solutior it this point must be close te itherw ise viditional 


sin must be used 
+ In reductions of pure nonradioactive D-glucom lactone Vields of approx 
mately 90 perce ver i | 
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Stress-Strain Relationships in Yarns Subjected to Rapid 
Impact Loading: 1. Equipment, Testing Procedure, 
and Typical Results ' 


Walter K. Stone, Herbert F. Schiefer, and George Fox ” 


Equipment is deseribed for elongating 
ranging from 10 to 100 meters per second 
100,000 to 500,000 precent per minute 


varns by 
The rate of straining at impact varies from about 


longitudinal impact at velocities 


\ procedure is discussed for obtaining load-elongation curves for loading and for unload- 
ing of the specimen and for loading to rupture in a time interval of only a few milliseconds 
The results of a typical loading and unloading test are presented 


1. Introduction 


Textile materials in processing and use are sub- 
jected ever more frequently to higher rates of 
straining. For example, in high-speed industrial 
sewing the thread is subjected to impact velocities 
ranging from | to 10 m/see as often as 5,000 times per 
minute, This results in high-frequency cyclic ac- 
celerations of the thread, which may equal several 
million centimeters per second per second, These 
high accelerations of the thread produce forces that 
approach, and may even exceed, the breaking 
strength of the thread. Similar conditions may be 
attained in high-speed processing of fibers, such as 
carding and combing or the weaving and knitting 
of varns into fabries; in the cord of an airplane tire 
during landing; in the safety line when a structural 
worker attached to it falls; and in the fabric, shroud 
lines, and webbing during the opening of a parachute 

The present work was undertaken to provide basic 
information on the behavior of textile materials 
when strained at very high rates; to relate the 
results to the molecular structure of the fibers and 
to the geometry of the varns and fabrics; and 
finally, to indicate the applications of this information 
to safety engineering, to industrial processing of 
textiles, to the development of fibers of improved 
impact properties, and to industrial and military end 
uses of textile structures for specific performance 
requirements 

This is the first of 
Impact testing In this paper 
procedures for this purpose are 
typical results are presented, 

The stress-strain behavior of textile and othe 
materials has been thoroughly investigated at rates 
of straining ordinarily referred to as “‘static’’ tensile 
testing. The velocity of the moving head of the 
testing machine in static testing is usually between 
about 0.1 and 10 inches per minute (4%107~° and 
t<107-* m/sec). The corresponding rates of straining 


au series of papers on rapid 
equipment and 
described and 


! rk por f t I Office of the Quartermaster General, Dey 


range from about 1 to 100 percent per minute for 
textile materials. In impact testing, usually made 
by dropping a mass a known height or by allowing 
a pendulum to swing freely through 90 to 180 degrees 
to rupture the specimen, the velocity at impact is 
usually limited by room dimensions to values between 
1,000 and 20,000 in./min. (0.4 and 8 m/sec). The 
corresponding rates of straining range from about 
10,000 to 200,000 percent per minute. Although 
excellent work has been reported in this velocity 
range {1 to 10],° the information concerning the 
stress-strain behavior of textile and other materials 
is not as complete as for static testing. The infor- 
mation at higher impact velocities, 20,000 to 
10,000,000 in./min (8 to 4,000 m/sec) or rates of 
straining ranging from about 200,000 to 100,000,000 
percent per minute, is incomplete. Much of the 
published work in this velocity range is theoretical 
or has been confined to very small strains [11 to 35). 
Meredith \7| has recently published excellent stress- 
strain relationships for viscose rayon, acetate, silk, 
and nylon yarns for rates of straining between 0.05 
and 70,000 percent per minute. The rates of strain- 
ing obtainable with the equipment described in this 
paper are between about 100,000 and 1,000,000 percent 
per minute for textile materials. Even this maxi- 
mum rate of straining is considerably below the 
rates of straining that prevail when textile materials 
are tested at their critical velocities. These rates 
of straining are estimated to range from about 
10,000 006 to LO0.000 000 percent per minute 


2. Equipment and Testing Procedure 


Kquipment has been constructed at the National 
Bureau of Standards for applying longitudinal im- 
pacts, al velocities up to 100 m/sec, to one end of a 
specimen of textile fiber, yarn, cord, or other ma- 
terials having similar dimensions. The rates of 
straining range from 100,000 to 500,000 percent per 
minute. A high-speed motion-picture camera and 
mirrors are used to photograph simultaneously the 





F igure n bracket nik the literature 
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positions of the head and tail ends of the specimen 
during the first few milliseconds from the time of 
impact. The stress-strain behavior of the specimen 
is then derived from the photographic record 

The specimen to be tested is usually 24 m long, but 
can be varied in length. One end ts fastened to a 
head mass, hereafter termed “head,” which is im- 
pacted The other end is fastened to a tail mass, 
hereafter termed “tail.” The head and tail, with a 
specimen of varn attached to them, are shown in 
figure | 

The manner of attaching the specimen is impor- 
tant. The specimen must not be damaged, high 
stress concentration must be avoided, and slippage 
must be minimized. An effective way of attaching 
& specimen is shown in figure 1. A small hole is 
drilled in the end of the head and the tail in such a 
way that it terminates in a helical groove on the 
surface. The bottom of the groove is coated with a 
rubber-like high-friction compound. The end of a 
specimen is threaded through the hole, wound in the 
and snubbed under the shoulder of a re- 
screw. It apparent that as the longi- 
tudinal tension in the specimen is increased during 
an impact test the radial pressure of the specimen 
against the groove also increased. The longi- 
tudinal stress is thus effectively transferred between 
the specimen and the head or tail through the in- 
creased contact friction. The effective length of the 
specimen is equal to the distance between the inner 
ends of the head and tail plus 1 em at each end for the 
helical groove 

The tail and attached specimen are inserted in a 
tube, as shown in figure 2. The right portion of this 
tube is made of transparent plastic, so that the tail 
and a portion of the specimen are visible. The head 
is fitted into the left end of the tube. Two disks are 


vroove, 
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Figure 1 Yarn specimen fastened to head and 


im pae i fest 


attached to the middle portion of the tube to serve 
as the bearing shaft when the tube is inserted in the 
impact equipment for a test. 

A portion of the impact equipment is shown in 
figure 3. The upper view shows the tube, with the 
head in a vertical position. Directly below the head 
and tube a portion of a heavy disk is visible. The 
disk can be rotated counterclockwise. Two lugs, 
called hammers, are attached to the disk. They are 
visible to the left of the head in the upper view of 
figure 3. In a test the speed of the disk may be 
adjusted to values between 5 and 50 rps. The 
peripheral speed of the hammers ranges from about 
10 to 100 m/sec. <A stroboscope is used to measure 
the speed of the disk. 

The upper view of figure 3 shows the hammers in a 
position where a sector (not shown) attached to the 
disk closes an electric contact in a triggering circuit 
When this circuit is closed a solenoid releases a spring- 
loaded cam. This cam rotates the tube containing 
the tail, specimen, and head through 90 degrees in 
less time than is required for one rotation of the 
hammers, even when the disk rotates at a speed as 
high as 50 rps. With the head in a horizontal posi- 


tion, lower view of figure 3, the head will be struck by 
the hammers of the disk for the impact test of the 
specimen attached to the head. 
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The assembly of the impact-testing equipment is 
shown in figure 4. A 5-hp variable-speed motor, 
used for rotating the disk at high speed, is shown at 
the leit Directly above it is a box in which the head 
is caught after impact. The tube containing the 
head, specimen, and tail is shown to the right of 
the box. Two seales, ruled in centimeters on trans- 
parent plastic, are shown directly above the tube. 
Photoflood bulbs are mounted directly in back of the 
A high-speed motion-picture camera is shown 
in the foreground. A system of mirrors is mounted 
the tube and camera. The paths of rep- 
resentative light beams from the head and tail are 
indicated by the broken lines. The mirrors are 
mounted on a metal base and can be readily ad- 
justed to different horizontal positions. The mirrors 
can also be adjusted in their inclination so that the 
image of the tail and its scale appear above the image 
of the head and its scale in the field of view of the 
camera. The head and tail are photographed in 
silhouette, as shown in figure 

After adjusting the speed of the disk to the de- 
sired value for a test, and with the head, specimen, 
and tail in position, the motor of the camera and 
the photoflood bulbs are turned on. When the film 
has attained the desired speed, the camera closes 
an electric contact in the triggering circuit of the 
solenoid This contact is in series with the solenoid, 
a series of batteries (alternating current of 60 cycles 
is not reliable for this purpose), and the contact on 
the disk referred to previously. The contact on 
the disk is closed the next time the hammers reach 
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Figure 5 Camera field view of specimen ends before im paet 
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the position shown in the upper view of figure 
At this time the solenoid is energized and actuates 
a trigger to release the spring-loaded cam. The 
cam rotates the tube and head into a horizontal 
position (bottom view of fig. 3). The head will be 
struck by the hammers on the next rotation of the 
disk for the impact test. The motions of the head 
and tail aiter impact are recorded by the high-speed 
camera, 

After impact, the head thrown to the left 
(fig. 4) and is caught in the box above the variable- 
speed motor. The tail is usually caught in the box. 
In some tests it remains in the tube or is thrown out 
of the rear of the plastic tube by the reaction from 
the specimen after rupture has occurred. 


Is 








3. Behavior of Specimen After Impact 


A few microseconds after impact, the head will 
attain a velocity exceeding that of the manner, and 
therefore will move in free flight ahead of the ham- 
mer. The sudden impact of the head produces a 
strain in the specimen where it is attached to the 
head. This strain depends on the impact velocity 
and on the material being tested. This strain and a 
tension approximately proportional to it will be 
propagated along the specimen toward the tail. 
When the strain reaches the tail it will be reflected 
and the reflected strain propagated toward the head. 
The strain and tension at the tail are approximately 
doubled upon reflection. 

The propagation of the strain along the specimen 
and the successive reflections at the head and the 
tail will continue. At each reflection the strain and 
the tension are increased. The increased tension 
will cause a deceleration of the head and an accelera- 
tion of the tail. This process will continue until 
the rupture strain is reached and the specimen 
breaks, or until the velocity of the tail equals that 
of the head, at which time the specimen will be at 
maximum average strain and under maximum aver- 
age tension, but both are below the rupture values 
After this time the tension in the specimen will 
continue to decelerate the head and to accelerate 
the tail, so that the velocity of the tail will exceed 
that of the head. As soon as the velocity of the tail 
exceeds that of the head the strain will be relieved 
and the tension in the specimen will be decreased. 
This of strain and propagation is 
probably similar to that cecurring during the in- 
crease of the strain, and will continue until all of 
the recoverable strain is relieved, at which time the 
tension will have been reduced to zero and the speci- 
men will become slack. The specimen is thus loaded 
to a maximum tension and extension and then un- 
loaded to zero tension and a set or extension in a 
total time of several milliseconds. Part of the set 
ma be reduced by delaved recovery. 

The propagation of the strain pulse after impact 
and subsequent reflections at the tail and head are 
indicated schematically in figure 6, A to H. Before 
impact, the specimen is represented in A by a series 
of adjoining circular springs. After impact, in B, 
the circular spring next to the head is deformed into 
an ellipse, and the strain pulse is indicated by the 
adjacent cross-hatched circle and ellipse. The 
arrow indicates the direction of strain propagation. 
C and D show the arrival and reflection of the strain 
pulse at the tail. At this moment the tail is set in 
motion, and the ellipse before reflection is deformed 
to one of greater eccentricity after reflection KE indi- 
cates the arrival of the reflected strain pulse at the 
head. After reflection at the head, the strain pulse 
is shown in F as traveling toward the tail about 
midway between the head and tail. In G the strain 
pulse is shown as traveling toward the head afte 
reflection at the tail and the high strain is indicated 
by the highly eccentric ellipses After the next 
reflection at the head it is assumed that the rupture 
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Schematic 
indicate ad by 
impact of the head at an instant between A and B 


FIGURE 6 re prese ntation of strain propagation an 


specimen, adjoining circular springs, afte 


Cross-hatched cirek the strain pulse 


ind ellipse represent 


strain is attained, and the ruptured specimen is 
shown in H. After rupture, the strain ts released. 
This condition ts represented in H by the 
version recovery of the ellipses of high eCCONLTICIES 
to ellipses of low eccentricity or to circles for com- 
plete recovery. The propagation of the strain- 
release pulse is indicated by the arrow. The 
ruptured end will snap toward the tail with a very 
high velocity, and the momentum acquired by the 
specimen, when the strain-release pulse reaches the 
tail, may be sufficient to stop the tail or to reverse 
its motion. 


recon- 


4. Tests on a Nylon Yarn 


A specimen of nylon varn was fastened at one end 
to a head weighing 19.82 g and at the other end to a 


tail of 4.85 g, the length between fastenings being 
65 em. The head was impacted at a velocity of 30 
msec, 


The motions of the head and tail were recorded 
photographically by means of high-speed motion 
pictures taken at a rate of 2,630 frames per second. 
The sequence of pictures is shown in ficure Z. be- 
ginning with two frames before impact, upper left- 
hand corner, and ending after release of the strain, 
lower right-hand corner. It can be seen that before 
impact in frames | and 2 the positions of the head 
and tail remain stationary relative to the fixed centi- 
meter scale, which is photographed on each frame 
After impact, which occurred between frames 2 and 
3, the head moves relative to the fixed seale until it 
finally disappears from the field of view of the 
camera. The tail, on the other hand, moves very 
little for several frames after impact. It then speeds 
up until in frame 10 it moves the same amount as 
the head. Thereafter it moves further per frame 
than the head; after frame 16 the distance moved 
per frame is maximum and essentially constant 
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Fiat RES Position of head and o} tail on successive frames 


in an impact test of a nylon yarn. 


The positions of the head and of the tail on each 
frame were carefully measured, relative to the fixed 
scale, with a precision micrometer. These positions 
are plotted in figure 8 against frame number and also 
against time in milliseconds. The resulting curves 
will hereafter referred to as distance-frame or 
distance-time curves. Zero time refers to the in- 
stant of impact and is obtained by extrapolating the 
distance-frame curve, A, of the head after impact 
back to the distance-frame curve of the head before 
impact. The point of intersection of two 
curves is taken as the instant of impact. 

The distance-time curve of the head is slightly 
concave downward, indicating that the head is slow- 


be 


these 


ing down. The slope of this curve at any point 
(time) is equal to the velocity of the head at that 
time. The velocity-time curve of the head is plotted 


in figure 9. Immediately after impact the velocity 
of the head is about 45 m/sec. At 4 msec after im- 
pact it is reduced to 30 m/sec. 

The distance-time curve of the tail, B of figure 8, is 
concave upward, indicating that its velocity increases 
with time. At 3 msee the slope of the distance- 
time curve for the tail is equal to that for the head, 
indicating that the velocities of the head and the 
tail are equal. The velocity-time curve of the tail 
is also plotted in figure 9. The velocity of the tail 
is zero before impact, then increases slowly immedi- 


ately after impact. It increases most rapidly at 3 
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FIGURE 9 Velocity of head and of tail after mpact in a test 


ofan jlon jarn 


msec after impact when it attains the head velocity 
of 35 m/sec. It then increases less rapidly becoming 
essentially a constant, 54 m/sec, at 5 msec. 

The slope of the velocity-time curve for the tail- 
figure 9, at any point (time) is equal to the accelera, 
tion of the tail. This acceleration is plotted against 
time in figure 10. The maximum value is more than 
a million centimeters per second per second. Be- 
cause the force on the tail (and thus on the specimen 
at the tail) at any time is equal to the mass times its 
acceleration, the curve in figure 10 is also a force- 
time curve of the specimen at the tail. The ordinate 
scale given at the right side of figure 10 is in force 
units (kilograms The curve in figure 10) thus 
indicates the force that the tail exerts on the speci- 
men during loading and unloading. The total time 
for loading and unloading is only 6 msec in this test 

From figure 9 it is seen that the deceleration of the 
head at 3 msee after impact is approximately equal 
to one-fourth of the acceleration of the tail. As the 
mass of the head is about four times that of the tail 
it is evident that the force at the head is essentially 
equal to that at the tail. 

In figure 8 the vertical distance between the dis 
tance-frame curve for the head, A, and the distance- 
frame curve for the tail, B, at any frame (time) is 
equal to the length of the specimen at that time 
The extension of the specimen at any time is given 
by the distance between the two curves at this time 
minus the distance between the two curves before 
impact. This extension is indicated graphically by 
the distance between curves A’ and B, where A’ is 
the distance-frame curve for the head, A, after it was 
shifted down to A’ so that the portion before impact 
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coincided with the similar portion of the curve B 
before impact. It can be seen that the extension 
(distance between curves A’ and B) increases for 
times up to 3 msec, and then decreases. It is zero 
again at 6.7 msec (frame 20). At this time, after 
impact, the head and the tail have moved the same 
total distance; that is, the length of the specimen 
(distance between head and tail) is again the same 
as before impact. If the specimen had recovered 
completely from its extension, it would have been 
under zero strain and zero tension at this instant 
It was observed, however, that the specimen was 
already slack in frame 19 (6.3 msec after impact 
The distance between the two curves at frame 19 
namely, at D, figure 8, therefore represents the ex- 
tension that remained after the specimen was com- 
pletely unloaded. Part of this extension may have 
been recovered subsequently by delaved recovers 
The extension during loading and unloading, ex- 
pressed as a percentage of the specimen length 
before impact, and termed ‘“elongation,”’ is plotted 
against time in figure 11. This elongation at any 
given time is the average elongation of the specimen 
Actually, the elongation is not constant along the 
specimen, but shows a sharp increment at the strain- 
wave front. 

From the force-time curve in figure 10 and the 
elongation-time curve in figure 11 it is possible to 
plot simultaneous values of force and elongation 
The result is a load-elongation curve for loading and 
also for unloading, as shown in figure 12. The 
energy to maximum load is then readily obtained, as 
well as the energy recovered during unloading. The 
difference between these values, ifconverted into heat 
and assumed to be uniformly distributed throughout 
the specimen, would be sufficient to raise the tempera- 
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lest on nylon yarn 

ture of the specimen about 10°C. Actually, this hys- 
teresis energy is not converted uniformly into heat 
throughout the specimen. For example, in another 
test in which the maximum strain was sufficient to 
rupture the specimen the fiber tips were fused, as 
shown in figure 13, which is evidence of high heat 
concentration at the point of rupture. This fused 
condition possibly results from the relatively low 
heat dissipation in nylon and the high localized 
energy released by bond rupture. 

The points for the loading curve, broken line in 
figure 12, appear to lie systematically above and 
below the solid straight line that is drawn through 
them. The straight line indicates that at this high 
impact velocity the nylon specimen essentially 
obeved Hooke’s law during loading. However, the 
scatter of the points about the straight line is too 
consistent to be ascribed to experimental errors. 
It appears that this scattering of the points about 
the straight line is experimental evidence that the 
force at the tail varies in a stepwise manner. Each 
step in load indicates the arrival and reflection of a 


strain wave at the tail. It is significant to note 
that the individual points of the velocity-time 


curve for the head, figure 9, show evidence that the 
deceleration of the head, and thus the force at the 
head, also varies in a stepwise manner. An analysis 
of these variations indicates that the time for the 
strain pulse to travel from the head to the tail, or 
0.65 m, is about 3107 see. The velocity of strain 
propagation is therefore about 2,200 m/sec. The 
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Ficgtre 12. Load at tail plotted against average elongation for 
loading and unloading of a nylon yarn in an impact test 








FIGURE 13 Magnified view of ruptured filaments of nylon 


show ng fuse d ¢ nds. 


strain at impact, given by the ratio of the velocity 
of the head at impact to the velocity of strain propa- 
gation, is thus about 2 percent. More detailed 
strain analyses for impact loading are given in 
subsequent papers [36, 37]. 
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Stress-Strain Relationships in Yarns Subjected to Rapid 
Impact Loading: 2. Breaking Velocities, Strain Ener- 
gies, and Theory Neglecting Wave Propagation’ 


Frank L. McCrackin, Herbert F. Schiefer, Jack C. Smith, and Walter K. Stone 


The behavior of a varn specimen fastened at one end to a head mass and at the other 
end to a small tail mass is analyzed for longitudinal impact of the specimen at the head. 
The analysis leads to a basic formula for “limiting breaking velocity,’ which is a character- 


istic property of the material and independent of the dimensions of the specimen. A simple 
procedure is described for obtaining its value. The values for cotton, nylon, and undrawn 
nvion varns tested at room temperature are found to be 130, 228, and 550 meters per second, 
respectively The practical application of the limiting breaking velocity to safety engineer- 
ing is shown by an example on safety lines. Formulas are derived for computing energy to 


any strain under impact loading, and the computed values are found to agree with those 
obtained from the area under the impact load-extension curves 


1. Introduction ; = 


The first paper of this series [1] * deseribes equip- 


ment and procedures for studying the stress-strain > a : J 
behavior of a varn specimen subjected to tensional 

impact In a typical test, a 65-cm specimen was Proure 1. Schematic drawing of test specimen 
attached at one end to a head mass and at the other or smpact tending 


end to a tail mass. The head was impacted at 
velocities between 10 and LOO m see, and subsequent 
behavior was recorded by high-speed photography 

The behavior of the yarn specimen after longi- 
tudinal Impact is treated theoretically in this paper, 
neglecting wave-propagation effects The results 


tually treated is that of a mass-less specimen to 
which a head mass (m+ 4)w and a tail mass (n+ 3)w 
are attached 


The force equations, neglecting air friction,’ are 


. » . Pr, 
are compared with those obtained experimentally for F=(n+4)w ye? 1) 
several materials Phe practical application of this 
work to safety engineering is indicated and 
| a 
fk (Wi 5 )W = 2) 
2. Theory dt 
Figure 1 represents a test specimen terminated by From definition of the center of mass, we obtain 
at head and a tail. The equations of motion of these n-+-4d)uus m+ 4 )wr,=(n+-m+1)wa 
masses are derived in terms of the following symbols Differentiating and eliminating w gives 
LL Length of specimen before impact 
Mass of specimet 1 dy 1 dy, | da ») 
? { : Mi > tT (Tit 5 | oe ot LL ee . >) 
Mass of tail a! dt os’ di dt 
\lass of head 
Position of tail before impact, relative to a fixed origin, O ; 
Position of head before Impact, relative to a fixed Because el. rs J (J's Dia), we obtain by 
_ origin, O differentiation 
lime after impact 
Position of tail at time, ¢, relative to a fixed origin, O j ] / 
Position of head at time, ¢, relative to a fixed origin, O ans ce —] naa (4) 
Position of center of mass of the svstem at time, /, dt dt “dt 
relative to a fixed origin, O 
«Strain in specimet » . 
/ Force acting on specimen at time, ¢ Solving eq (9 and (4) for di dt and dy; dt gives 


Velocity of head immediately after impact 


dr, da m+43)L de 
It is assumed here that aw is concentrated at dt dt , } + dt’ ae 


the head and the remaining 3 is concentrated at 


the tail of the specimen. Thus the problem ac- and 
~nlanininai dr, dz., (nt+43)L de ' 
Nork V ponsore the flice of the uarterm ter reneral epart T ° ») 
the aon y, and pr ciao I he at ‘the Septen Ne 1953 he dt dt nm) dt 
Che Fiber Societ nd a it the November 1954 meeting of the 
Rheology Che force on the tail due to air friction was calculated to be approximately 
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At maximum strain e=e,, and de/dt=0, therefore, 
dr, dt=dr, dt=dr,/dt. That is, for a typical loading 
and unloading test under impact conditions, the 
velocity of the tail mass at the instant of maximum 
strain is equal to the velocity of the head, and both 
are equal to the velocity of the center of the total 
mass. <As the total momentum imparted by the 
hammer to the head immediately after impact is 
m-+4)wr, it follows from conservation of momen- 
tum that = dr,/dt=(m+})ry/(n+m+1 constant, 
which is also the velocity of the head and of the 
tail at maximum strain. The velocity of the head, 
therefore, decreases from ¢) immediately after im- 
pact to (m+4)r/(n+m+1) at the time of maxi- 
mum strain. It is seen that this decrease is equal 
to n- , Yo/(n+m+l and = thus depends upon 
the values of n and m 

From eq (5) and (1), or (6) and (2), it follows that 


t- n+4)(m+ 3 Lu 


Integration gives 


dle » jes" he | "e : ’ 
( = ) Lurin- : m- | Fide ré 
At f (), € 0, and de dt ve/L Therefore. C 
L and 
ile : / 2 nm) Pies 
( dt ) L ; Lu lin 4 jin } | Fide 10 


At the maximum strain ee, and dedt—0 lt 
follows. therefore. that 


1 1 é *. 
Vi > Mint» oy 


log W log _ = +log "Fide 1] 
~ nami Vw. 


In an impact test the specimen will break if the 
impact velocity, 7, is large enough that the rupture 
strain, e,, is attained. The minimum breaking ve- 
locity, ¢,, for a given n is the velocity just sufficient 
to cause the specimen to break when ¢ é,. Under 
these conditions we have 


vi -7 } in , 2L ‘ , 


tlogy =| File 12 


Equation (12) expresses the minimum breaking 
velocity as a function of n If «, and the shape of the 
stress-strain curve are both independent of the rate 
of straining, the expression under the radical sign is 
a constant having the dimensions of velocity squared 
At high rates of straining the value of ¢, is less than 
that obtained at conventional testing rates. Also, 
the slope of the stress-strain curve is greater at higher 
rates of straining. However, for rates of straining 
of the same order of magnitude, these changes are 
small, and their effects compensate each other C'on- 
sequently, the area under the stress-strain curve 
should be essentially constant for the rates of strain- 
ing considered here. Under these conditions a plot 


of log rv, versus log(n+4)(m+4) (n+m—+1) would 


log r t log 





be a straight line of slope— } 
By extrapolating from impact test data to the 


point at which log(n+4)(m+3)(n+m+1)=0, a 
characteristic velocity, r,, is obtained: namely, 
2] *. 
= File 1:3 
\ ©. v 


If we let p densit of the specimen and o stress, 
eq (13) reduces to the form 


5 / 

\ = ale 14 
PJo 

It is thus seen that 7, is a quantity characteristic of 

the material itself, except for a possible dependence 

on the rate of straining 


Equation (14) can be rearranged in the form 


pr ade, 1D) 


which states that the kinetic-energy density in the 
specimen when traveling at velocity ‘ Is just equal 
to the strain-energy density required to break the 
specimen. We thus see that if the specimen is 
impacted at a velocity greater then 7,, it will always 
be broken. For this reason, we call v, the limiting 
breaking velocity 

Von Karman [2] derived a critical velocity, 7, at 
which a filament will break immediately upon 
Impact 

me | do 


J. Wea” " 


This equation Was obtained from a consideration ol 
plastic and elastic wave propagation in a material 
for which the stress-strain curve is concave down- 
ward, and the stress ts independent of the rate of 
straining 

The results calculated from eq 14) and (16 
differ slightly. However, in the special case when 
Hooke’'s law is obeved both eg 14) and (16) reduce 


to 


yp’ =e li 


This formula may be expressed in terms of the 
more familiar textile quantities of tenacity and 
percent elongation at rupture, namely, 

vy’ 29.74 tenacity, elongation, 1s 
where 7’ is In meters pel second, tenacity is in grams 
per denier, and elongation is in percent 

Meredith [3] has obtained evidence that at high 
rates of straining the load-extension curves of som 
materials become more Hookean. These formulas 
thus assume considerable lmMportance and by means 
of eq IS estimates of the eritical velocities ol 
varns may be computed from their tenacities and 
elongations at rupture under lmMpact testing condl- 
Lions 
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In a subsequent paper [4] the behavior of the 


; 
‘ 


specimen after longitudinal impact is treated by a 
theory that the effects of wave propaga- 
tion. This treatment shows that the results given 
in this paper for nonwave theory are valid when the 
tail mass is greater than 40w. For tail masses less 
than 40w, the effects of wave propagation may 
become appreciable and must be considered in’ the 
theoretical treatment 


considers 


3. Experimental Results 


experimental values of 1 corresponding to differ 
ent values of nm were obtained on a number of dif 
ferent materials These values of log ry, are plotted 
in figures 2, 3% $ against ; +) 
7, ]) These plots are straight lines having 
the empirical equations 


ana log (n m 


Mi 


! " 
5 ‘it 
log ) 0.496 log = “’--Llog 130 
4 ‘it 
for cotton thread 
= Vi i Mt , 
low ) (0.495 log . « og 228 
Vi ‘ti | 
for nvlon varn and threads 
1 1 
Vi ; ti > — 
log ; 0.500 log = +Jlog pou 
fi ‘ti | 
for undrawn nvlon 
It should be noted that the slope of each line 
represented by thes equations is very close lo 


in accordance with eq (12) for log ; 


4. Design of Safety Lines 


The empirical equation derived for log / has many 
practical applications. Consider, fos example 
safety line of length LL. that connects a construction 
worker to a mgid point of the Under 


these Infinite The maximum free 


structure 


conditions, m 


is 


279 


400 











@ ONE PLY NYLON 
t QO TEN PLY NYLON 
V NYLON THREAD 200 ¢g 
) BONDED NYLON 4 
‘ A SONDED NYLON t 
2 00 & 
* Bo > 
o , 
QO % VU 
4 iss 50 Go 
n 
. wi 
a! 
bad Oy 
va) (r+! . é 
LOG V.+-—0.495 LOG s+ ims 72) tLOG 228 . 
n nemei 
! a ee ee ee ee Pcustbanaticantivendinndl 4 id 
0 i 2 
Loe (nt 4)im+) 
nemei 
bict Kt 3 Relationsl p hetween impact velo a iplu 
i and head-tail mass paramete " ; ” ; ’ ” / 
for nulon yarns and threads 
5 
UNORAWN NYLON 
1 
Vavir v7 
LOG V,*- 0.50 LOG eeu & +106 550 ; 
| 200 & 
e 
c 2 100 > 
£2 
> 
> ; 
| , 
1 $0 ° 
| ad 
5 | 
1 20 
| 
i sacihslialacdiihieailesiiae oe 
ine Ze)im+/e) 
GG Eee . 
neme 
hrourt } Relationship heltween smpact velocity to ruplu 
i and head-tail mass parameté 7 5 m a) /An m / 
for undrawn nylon 
fall of the worker would be 24, and the maximum 
free-fall velocity 2 WEP As m becomes infinite 
eq (12 reduces lo 
1 
log 7 t log(n + 4) +log 7 14 


be 


behaves 


\Iaking the ussulmption (which must checked by 
experiment that the safety line unde 
pact like the single yarn considered in the deriva 
19 


one obtaims 


tion of substituting + 2,91 and solving 


for 


eg 
fi 
) 


=] 
ty 


L 


Mi 


If vy, is 200 m see for the material of the safety ling 
L—10m, and g=9.8 m sec*, the computed value of x 
is 101. If the weight of the worker, nw, is 202 Ib, the 
minimum weight of 10 m of safety line, w, to just 


support the worker at is 2 lb 
this computation 

This weight is approximately that of a rope 
in diametes A nylon rope i) 


breaking strength of 4,000 Ib 


Impact according to 
It) 
in diameter has 
a) Jt) Linnie. 


5 t 
stutie 


il ol 











the weight of the worker. However, under the 
dynamic conditions just assumed, a nylon rope with 

200 m/sec would have no factor of safety what- 
ever according to this computation. 

If », for the material of the safety line ts 100, 
n—=25, and for the same weight of worker, w would 
have to be 8S lb. This weight is approximately that 
of a rope about 1 in. in diameter [5]. A manila rope 
| in. in diameter has a static breaking strength of 
11,000 Ib, or 50 times the estimated weight of the 
worker. Again, under dynamic conditions, this 
manila rope would not provide a reasonable factor of 
safety according to this computation. Results 
obtained by Newman [6] showed that a %,-in.-diam 
sisal rope 10 ft long broke when a weight of 142 Ib 
was dropped a height of only 20 ft. 


5. Energy of Straining 


The energy per unit length, (?, to stress the speci- 
men to a strain eis given from eq (10) by the equa- 


tion 


*e ' Law (n , m , , le ‘ 
W hile = / ee ( r)—( a) 2() 


where ded? is the slope of the strain-time curve at 
If the specimen is not ruptured but 


, and then unloaded, 


the strain € 
is loaded to a maximum strain, ¢€ 


the energy, (,,, to the maximum strain, «,,, IS given 
bv the equation 
e, ] 1 
Lu Vi 5 ‘ti 3 D 
() hile = ( ). Z 
. , 2 di Mii | L 
since (de dt is zero. If the specimen ts ruptured, 


the energy, (,, to the rupture strain, €,, is given by 


the equation 


" ' Lu (n ; nl 4 , ile 
( File = | : ( L ) —( dt ) 


Hi Mi 


The value of (de/dt),., will be zero, if the impact 
velocity of the head is just sufficient to produce 
rupture. Above this velocity (de/dt),.. is a posi- 
tive quantity 

By means of eq (20), (21), and (22), it is possible 
to compute (, &/ and (). from the conditions of a 
test For it Is only necessary to know the 
velocity, “, of the head at impact. For @ and 
the value of de/dt and (de/dt),.. must be known in 
addition to the value of ¢ The values of 7, ded. 
and (ded), are readily obtained from the posi- 
tions of the head and of the tail on the high-speed 
motion pietures taken during the test The energy 
values computed by means of these equations for 
different materials and impact testing conditions are 
plotted in figure 5 against those derived from the 
area of the impact-load-extension curves. The 
plotted points fall near a straight line having a slope 
equal to |. 
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6. Summary and Conclusions 


The theoretical treatment given in this paper for 
longitudinal impact of a varn leads to a formula for 
the limiting breaking velocity that is a characteristic 
property of the material and is independent of the 
dimensions of the specimen. A specimen impacted 
at a velocity greater than its limiting breaking 
velocity will always be broken. A simple procedure 
is deseribed for obtaining this limiting breaking 
velocity from impact-test data, and results are given 
for cotton, nvion, and undrawn nvlon Varns The 
potential importance of the limiting breaking velocity 
to safety engineering is indicated by its application 
to estimates of margin of safety of rupture of safety 
lines for censtruction workers The theoretical 
treatment also leads to formulas for obtaining 
energy to any given strain, including rupture strain, 
under different impact conditions. These values 
agree reasonably well with those obtained from the 
area under the impact-load-extension curves 
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Waterdrop Collisions With Solid Surfaces 
Olive G. Engel 


An approach is described to the difficult problem of the dynamies of 
The flow in the impact plane is traced chemically. 


sphe re 


ill impinging water 
The stages In + collapse of the 


water mass and of its radial flow are shown in photographs that were taken with a high-speed 


motion-picture camera 
force and of the radial flow 
the radial flow is considered 


velocity 


of the water after the collision. 


l. Introduction 


Erosion by waterdrop impact has recently received 
attention because of the damage that is done to high- 
speed aircraft on flying through rain. The problem 
is not new. Similar erosion is produced in steam 
turbines by the waterdrops in wet steam. 

The destructive force causing this erosion results 
from the collision of the solid surface with the water- 
drop. At high impact velocities a waterdrop acts 
as though it were a hard sphere. Unlike a sphere 
of hard material, it undergoes an outward radial flow 
of very high velocity as a result of the collision. 

The question as to the order of magnitude of the 
pressure that is developed in collisions of waterdrops 
against solid surfaces has been discussed since the 
first studies of this erosion were begun approximately 
30 vears ago. Efforts have been made both to ealeu- 
late the pressure from theoretical considerations and 
to measure it experimentally. 


2. Behavior of a Liquid Drop on Collision 
With a Flat Solid Surface 


The process of collision of one solid against another 
solid is well known in elasticity theory. The collision 
of small solids with a body of liquid has also been 
investigated in water-entry problems. The collision 
of liquid drops against a solid surface, however, has 
been studied very little. Some preliminary work for 
the purpose of obtaining a better understanding of 
waterdrop-to-solid collisions is described in- this 
paper. This work includes (1) a chemical mapping 
of the radial water flow in the impact plane, (2) the 
use of high-speed motion-picture photography to 


“stop” the motion of the drop sufficiently during thi 
collision so that the stages in the transition from 
vertical to radial flow can be observed. and (3) the 


use of schlieren photography to study details in the 
radial flow. The following discussion applies to the 
case Of collision of a waterdrop with a flat solid 
surface. 


2.1. Map of the Radial Water Flow 


Worthington [1]! with the aid of spark photogra- 
phi made a study of the forms assumed by drops of 


I irt n brackets indie iture references at the end of this paper 


Empirical determinations of the 
are reported The 
A semiempirical 
assumptions leads to equations for the maximum impact pressure 


time dependence of the 
possible occurrence of cavitation in 
analysis based upon various simplifying 
and for the rate of spread 


Impact 


liquids (mainly milk and mercury) after a vertical 
fall from heights up to 11 in. onto a horizontal sur- 
face. A tendency toward formation of radial arms 
was found by Worthington to increase both with the 
height of fall and with the size of the drop. 

The radial flow of a waterdrop was mapped chem- 
ically at the National Bureau of Standards. To 
accomplish this, a very small crystal of sodium di- 
chromate held on the point of a needle was placed in 
the bottom surface of a waterdrop just before it fell 
from a flat-nosed pipet. The solution of this oxidiz- 
‘ng agent is heavier than water and remained in 7 
bottom of the drop. The drop was then allowed t 
fall onto a glass plate covered with a filter paper that 
was previously wet with acidified starch and potas- 
sium iodide solutions. A typical starch-iodide print 
of the radial flow from a waterdrop that fell approxi- 
mately 1.5 ft is shown in figure 1. The water that 
struck first, and which contained the sodium dichro- 
mate, washed to the periphery of the flow, as can 
be seen by the stronger starch-iodide color there. 
The water that struck last essentially did not flow. 





Starch-iodide trace of the 
after collision with a solid surface 


FIGURE 1 radial flow of a wate rdrop 
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The radial arms in the structure show that the flow 
was channeled. It is conceivable that these channels 
are a result of the viscous drag between the solid 
surface and the flowing liquid, as was suggested by 
Worthington [2]. On the other hand, the observa- 
tion of Worthington that the number of channels 
increases both with the size and with the height of 
fall of the drop suggests that the number of them 
may be related to the kinetic energy of the drop 


2.2. Stages in the Collapse of a Waterdrop After 
Collision 


High-speed motion pictures can stop sufficiently 
the movement of a waterdrop after it has collided 
with a solid surface so that the horizontal flow can 
be observed. Figures 2, 3, 6, 7, and 9 are enlarge- 
ments from high-speed moving pictures that were 
taken by T. C. Hellmers, Jr., and Harold Goodman, 
both of the Diamond Ordnance Fuze Laboratory, 
and by Warren P. Richardson of the NBS Technical 
Reports Section. The camera used is capable of 
taking 15,000 frames per second. The operating 
velocity at which any incident is photographed can 
be computed from the 60-cycle timer marks along 
the border of the film. 

The steps in the collision and flow of a waterdrop 
after free fall through distances of about 1.5 ft and 
of about 20 ft are shown in figures 2 and 3, respec- 
tively. For the waterdrop collisions shown in these 
figures a very small crystal of sodium dichromate 
was placed in the bottom surface of the drop, and 
the drop was allowed to fall against a glass plate 
covered with a filter paper that was previously wet 
with acidified starch and potassium iodide solutions. 
A waterdrop oscillates between a vertical and a 
horizontal elongation as it falls. The collision 
shown in figure 2 occurred while the drop was elon- 
gated vertically; that shown in figure 3 occurred 
while the drop was elongated horizontally. It can 
be concluded from these pictures that the bead of 
liquid of the drop, before it enters into radial flow, 
shows an ability to resist a change of shape during 
collision with a solid. This may be a result of its 
inertia, or of its viscosity, or surface tension. 


2.3. Schlieren Pictures of Impinging Waterdrops 


To accelerate waterdrops to their terminal velocity 
in air by free fall under the force of gravity, a tube 
was installed in a staircase well to extend through 
three stories of the Industrial Building of the Na- 
tional Bureau of Standards. This tube served the 
purpose of shielding the waterdrops from air cur- 
rents so that they could be made to fall with little 
deviation on a predetermined spot on a glass plate 
located at the lower exit of the tube. <A flat-nosed 
pipet, from which the drops originated, was mounted 
at the top of the tube. 

It the solid surface against which the waterdrop is Figure 2. Stages in the flow of a waterdrop that fell from a 





allowed to impinge is a glass plate, the stages in the height of approrimately 1.5 fee 
radial flow of the drop can be photographed through 
the plate itself. A complete schlieren svstem was 


found to produce pictures from which the most infor- 
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FIGURE 4 Elements of a schlieren system for photographing 
impinging waterdrops 


mation could be obtained. A sketch of the elements 
of the arrangement that was used is shown in figure 
$4. The arrangement was planned by Roland V. 


Taek ore Shack of the NBS Optical Instruments Section. 


a. Schlieren Arrangement for Photographing an Impinging 
Waterdrop 


The light source was a carbon are searchlight that 
operated on a 40-amp current. The light was colli- 
mated by a parabolic mirror and was directed down 
the tube, through which the waterdrops fell, by 
means of a front-surface plane mirror. The tube 
was dved black on the inside to reduce light seatter- 
ing. The degree of collimation of the light was 
further improved by a circle of cardboard that was 
fastened to the exit of the tube and that allowed only 
the central core of the collimated beam to pass 
through to the glass plate. 

The collimated light was transmitted through the 
glass plate and was focused by a lens. The sharp 
spot of light, which formed at the focal point of the 
lens, was received on an opaque barrier that had the 
same size and shape as the spot of light itself. This 
barrier was made by allowing the spot of light to 
fall on a spectroscopic plate that was rigidly mounted 
at the foeal point of the lens in complete darkness. 
The opaque barrier served the purpose of a schlieren 
knife edge. The camera was mounted below. the 
opaque barrier and was focused on the upper surface 
of the glass plate. Theoretically, only the light that 
was scattered by the impinging waterdrop would be 
deflected so as to enter the lens of the camera. 
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b. Interpretation of the Schlieren Patterns 


a To show how the light that was scattered by an 
height of approximately 20 feet impinging waterdrop would appear in pictures taken 
with the schlieren arrangement, four glass models of 
waterdrops were fashioned by L. Testa of the 
NBs Glassb lowing Shop. The models ranged 


, 


j 


736155 283 








from a solid glass sphere to an almost flattened glass 
disk. They are shown in figure 5, along with the 
corresponding photographs that were made when 
they were placed on the glass plate of the schlieren 
arrangement. The pictures numbered 1 through 
4 in figure 5 are the glass models. The pictures 
having the corresponding primed numbers are the 
schlieren patterns that the models produced. <A 
50-mm lens was used to photograph the schlieren 
patterns. The lens was kept at its widest aperture 
so that the depth of field was very short. 

In picture 1 of figure 5 a hollow glass bead (black 
pointer) is shown with the solid glass bead. The 
schlieren patterns obtained from these beads show 
how difficult it is to tell the difference between a 
bubble in water and a protuberance from the water 
surface. The hollow bead produces only a somewhat 
larger high light in the schlieren pattern than the 
solid bead produces. Picture 2 of figure 5 shows a 
glass model of a waterdrop in the first stage of radial 
flow. The schlieren pattern produced by it is the 
same as that produced by the solid glass sphere but 
has a slightly larger diameter. A small difference in 
diameter is meaningful in this case because the 
geometry of the schlieren arrangement remained 
constant. 

Pictures 3 and 4 of figure 5 show glass models of a 
waterdrop in more advanced stages of radial flow. 
The glass model in picture 4 has a ridge at the periph- 
ery of the radial flow, whereas that shown in picture 


3 does not. Comparison of the structure of these 


models with their sehlieren patterns shows that a 
flat area appears light in the schlieren pattern, regard- 
less of whether it exists on a ridge or in a depression, 
whereas an area that has anv degree of inclination to 
the reference plate on which the 
appears 


model is 
the only way 


glass 


resting dark. Apparently 





models ot waterd Ops 


Glass 


Figure 5 


and the 


that the presence of a ridge or a depression can be 
identified is by the sloping area leading to it. The 
spherical residue in the center of the glass models in 
pictures 3 and 4 represents an inclined area and 
appears dark in pictures 3’ and 4’. In pictures 1’ 
and 2’ the schlieren patterns are also dark as is now 
expected because the glass models that produced 
them (pictures 1 and 2) have no area that does not 
have an inclination to the reference plane. 


c. Increase in Size of the Central High Light 

The schlieren patterns shown in pictures 1’ and 2’ 
of figure 5 have a central high light as do also those 
shown in pictures 3’ and 4’. This results because 
the glass sphere or spherical residue acts as a lens. 
If the spherical residue should become progressively 
more flattened, it should behave as a lens of longer 
focal length, and the high light should be seen to 
become larger and more diffuse. Although this 
effect is observed in schlieren pictures of the collision 
of real waterdrops with the glass plate, it is not 
observed in pictures 1’ through 4’ probably because 
the degree of curvature of the spherical residue in 
the glass models is not sufficiently different. 

The first 10 frames of a collision incident of a real 
waterdrop impinging against a glass plate as photo- 
graphed at a small oblique angle above the glass 
plate are shown in figure 6. As soon as the water- 
drop contacts the glass and flows at all, it becomes 
a combination lens in which the lower lens has a 
small central flat area. Under this condition two 
foci exist, namely, that of the double-convex lens 
and that of the single planoconvex lens. <As_ the 
solid surface moves further through the waterdrop 
the intensity at the focus for the double-convex lens 
decreases and that at the focus for the single plano- 


convex lens increases. When the head of water that 





primed n imbers 


schlieren patterns 
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emains of the drop has decreased to a hemisphere, 
mly the focus for the single planoconvex lens exists. 
The focal length of a sphere of water measured 
from the point of contact with the flat surface on 
vhich it is resting, for rays making a small angle with 
the axis, is equal to its radius. The radius of curva- 
ture of the top of the drop in the first frame of figure 
6, as found by trial and error with dividers, was 
0.36 em. The focal length of a single planoconvex 
lens of water measured from the plane side, which is 
in contact with the flat surface on which it is resting, 
for light that is incident on the curved surface of it, 
is 2.25 times its radius of curvature. As the radius 
of curvature of the water hemisphere that remains of 
the drop in the fifth frame of figure 6 was 0.51 em, 
its focal length at this stage of the collision was about 
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1.15 em. The difference in the focal length of the 
water lens as it exists in frame 1 and in frame 5 is 
0.79cem. As the depth of field of the camera used to 
take the schlieren photographs was less than 0.2 cm, 
it seems reasonable that a perceptible increase in the 
size of the high light in the center of the schlieren 
patterns of a waterdrop should be observed during 
the course of a collision between the waterdrop and a 
glass plate. This effect can, indeed, be observed in 
frames 1 to 13 of figure 7. 


d. Lifetime of the Head of Water of a Drop That Collides at Its 
Terminal Velocity in Air With a Glass Plate 


The glass models were placed directly in the center 
of the schlieren arrangement in line with the opaque 
barrier and with the center of the camera lens. The 
waterdrops were usually displaced from this position. 
This may explain why the schlieren patterns of the 
waterdrops show an additional high light on the side 
of the head of water of the drop as long as it exists as 
such. When the head of water just disappears so 
that only a disk of water in radial flow remains, the 
side high light and the central high light merge. For 
waterdrops impinging against a glass plate at their 
terminal velocity in air, the number of frames from 
the first point of impact to the merging of the two 
high lights indicates that the time required for the 
disappearance of the head of water of the drop is 
about 1 msec (millisecond). 

The merging of the high lights seems to be accom- 
panied by a flash of light, which unfortunately ob- 
scures the details in the radial flow for a space of 6 or 
7 frames. The explanation of the light flash seems to 
be that the water lens at this stage of the flow of the 
drop acts in conjunction with the glass lens of the 
schlieren system to allow light to pass around the 
opaque barrier. The event occurs consistently. 


e. Structures in the Radial Flow 


On inspection of schlieren patterns produced by 
collision incidents between a glass plate and imping- 
ing drops, such as those shown in figure 7, it is 
seen that the radial flow of the liquid appears dark. 
With reference to the observations made in regard to 
the schlieren patterns produced by the glass models, 
this appears to indicate that the radial flow of the 
liquid has some degree of inclination with respect 
to the glass plate. Worthington [1] made observa- 
tions with drops of milk and mercury for low heights 
of fall that show that the liquid in radial flow is 
deepest at the periphery (see fig. 15). 

In the schlieren pictures the periphery of the 
radial flow is bounded by a beady white line. The 
beady character of the periphery of flow can be 
detected as early as frame 5 in figure 7. With 
reference to the schlieren patterns produced by the 
glass models, a white area is either the top of an 
elevation or the bottom of a depression. In the 
case of the white boundary of the radial flow it is 
most likely an elevated ridge. This conclusion 
seems to be supported in the pictures shown in 
figure 2. 
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FIGURE 8S 


Ridges formed in the flow of a de forming lead pe let 


atler collision 


Waves that run out over the radial flow can also 
be seen when high-speed pictures of a collision of a 
waterdrop with a glass plate are viewed with a 
projector, This effect is not visible In still 
enlargements 

Finally, inspection of pictures of various collision 
incidents between waterdrops and a plate 
shows that structures of some kind develop in the 
radial flow. The development of these structures 
follows the same general course. The central high 
light increases in size and finally takes on a definite 
structure. The structure grows in complexity and, 
after the central and side high lights of the drop 
merge, finally appears as separate white spots. That 
these white spots or areas are probably ridges and 
depressions can be seen from the flow that a de- 
forming lead pellet undergoes after impact at very 
high speed against a metal plate. The flow of a 
deforming lead pellet is shown in figure 8. In the 
case Of the radial flow of water after impact, these 
ridges eventually vanish. 

The explanation of why these ridges should form 
is not clear. The compressional wave started by the 
Impact moves through the drop at the speed of 
sound in water and reflects from the opposite side 
of it as a tension wave. Thereafter it changes 
from tension to compression and back to tension 
alternately at the end of each round trip through the 
head of water that remains of the drop. It is possible 
that a standing wave may be produced as a conse- 
quence of the boundary conditions imposed on this 
alternate wave of compression and tension. On the 
other hand, it is possible that a solution that would 


glass 


predict standing waves may result from the con- 
tinuity equation and the boundary conditions of the 
radial flow in the impact plane. 


f. Cavitation in the Radial Flow 


One of the most difficult aspects to understand in 
the problem of the erosion of durable structural 
materials as a result of high-speed waterdrop 
impingement is how initiation centers develop. It 
has been suggested by Albrecht Herzog of the Wright 
Air Development Center that incipient cavitation 
erosion may provide a first surface roughness from 
which erosion that results from waterdrop impinge- 
ment could progress. <A first surface roughening is 
observed both in the case of cavitation erosion and 
in the case of erosion by waterdrop impact. From 
this standpoint it is of interest to know whether or 
not any evidence of cavitation can be detected in the 
radial flow of a waterdrop after collision with a solid 
surface. 

Cavitation is the formation of bubbles in a liquid. 
[t occurs When the pressure on a liquid or in a small 
volume in a liquid drops below the value of the vapor 
pressure of the liquid at the temperature in question, 
The cavitation itself is the spontaneous 
srowth of gas nuclei (which already exist in the 
liquid) due to vaporization of the liquid across the 
liquid-vapor interface that each growing nucleus 
provides. When the pressure on the liquid is raised, 
or when the bubbles move out of a local low-pressure 
region in the liquid to a region of high pressure, the 
bubble-cavities collapse. It is the collapse of the 
cavities that produces the type of damage known as 
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cavitation erosion, 

There are at least two wavs by means of which it is 
possible to explain theoretically how the conditions that 
produce cavitation may develop in a liquid drop after 
it collides with a solid surface at a velocity of suffi- 
cient magnitude. The first possible explanation is 
based on the very rapid radial flow of the liquid 
It suggests that when the head of liquid of the drop 
has just disappeared into radial flow, the continued 
outward flow of the liquid under its own monemtum 
will produce a drop in pressure at the center of the 
spreading liquid disk. If the pressure at the center 
of the spreading liquid disk should fall below the 
vapor pressure of the liquid in question, cavitation 
may occur. This explanation has been advanced by 
Dr. Herzog. That a liquid drop may break in the 
center of the radial flow is shown in the work of 
Worthington [1]. He found that impinging drops 
of mercury break in the center and form a ring rather 
than a solid disk of liquid. 

The second possible explanation of how cavitation 
conditions may be realized is based on the alternating 
wave of compression and tension that exists in the 
head of liquid that has not vet become part of the 
radial flow of the drop. At the first instant of col- 
lision a wave of compression, which is initiated by the 
impact, moves through this head of liquid to the top 
of the drop. Here it reflects from the free liquid-to- 
air surface as a tension or negative pressure wave. 
This negative pressure wave moves back through the 
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drop to the impact surface. It must, in fact, be 
focused to a very small area of the impact surface by 
the curvature of the top of the drop in much the same 
way as ultrasonic waves in water may be focused by 
a watch glass. The returning negative pressure wave 
adds algebraically to the compressional wave that is 
still being initiated at the impact surtace as a result of 
the collision. Since, however, the impact pressure 
has been decreasing steadily during the time interval 
in which the first compressional wave front made its 
round trip through the waterdrop and returned as a 
negative pressure wave, the net pressure is negative. 
The extent to which the net pressure is negative 
depends on the rate of decrease of the impact pressure 
with time. This explanation of how cavitation con- 
ditions may be produced in a liquid drop after 
collision with a solid surface is due to the late Francis 
EK. Fox of the Catholic University. 

Some slight evidence was found for what might be 
cavitation in schlieren pictures of a waterdrop col- 
lision with a glass plate. At the suggestion of 
Virginia Griffing of the Catholic University, water- 
drops saturated with argon gas were used in an effort 
to show that the evidence that was observed was 
actually produced by bubbles in the water. Dr. 
Griffing has found that cavitation occurs readily in 
water that has been saturated with argon gas 

Enlarged pictures of 75 consecutive frames from 
one of the collisions of a drop of argon saturated 
water with the glass impact plate are shown in figure 
9. In this figure the life evele of a bubble can be 
followed from frame 19 to frame 33. The appear- 
ances, Which are described here, are more readily seen 
in the original photographs than in the halftone re- 
productions. The bubble appears slightly to the 
right of and below the center of the radial flow in 
frame 19. It grows in size and distinctness, then 
diminishes, and eventually vanishes in frame 33. As 
the time per frame for this collision incident was 93 
usec, the lifetime of this bubble was 1.4 msec. This 
is a reasonable lifetime for a cavitation bubble [3]. 
This bubble made its appearance shortly after the 
head of water of the drop vanished into radial flow in 
frame 14. 

A second bubble life cvcle can also be followed in 
figure 9. This bubble appears in frame 24 on a level 
with the bubble that has just been discussed but 
about ‘s in. from the periphery at the left side of 
the radial flow. It increases in brightness and then 
dims virtually to extinction in frame 34. The life- 
time is 0.9 msec. In frame 42 a white spot, which 
may be the same bubble, starts to increase in bright- 
hess. 
the radial flow. It continues to increase in bright- 
ness and, at frame 60, that is, 1.7 msec aiter its 
reappearance, it drifts into the left periphery of the 
radial flow. 

Two bubbles slightly below and to the leit of the 
center of the radial flow can be observed from frame 
31 to frame 45. In frame 46 the characteristic white 
semicircle which seems to be associated with a surface 
bubble (possibly due to the raised rim of water that 
accompanies such a bubble) appears. One of the 


It is about \%. in. from the left periphery of 


‘ :; 
bubbles has@broken by frame 52, as is seen by the 
train of ripples that forms from the white semicircle. 
The other bubble becomes obscure at about frame 


59 or 60. A bubble floating on a liquid surface de- 
presses the liquid level below it and is accompanied 
by a stable meniscus or raised rim of liquid at its 
periphery. If the bubble collapses, the differences 
in liquid level are no longer stable. This condition 
will produce waves. These waves are predominantly 
capillary waves since the wave length is small. 

The collision incident for a drop of ordinary hy- 
drant water against the glass impact plate is shown 
in figure 7. A small white spot can be seen very 
slightly to the right of the center of the radial flow 
inframe 14. This spot grows in brightness to about 
frame 25. It then becomes progressively dimmer 
until it disappears at about frame 40. This 
lifetime of about 2.5 msec. The rise of a bubble to 
the surface and its collapse in the surface can also 
be traced in this figure. The first evidence of this 
bubble is in frame 29 at about the center of the 
radial flow. By frame 42 it has developed into two 
white semicircles enclosing a small white spot. As 
has already been noted, this structure seems to be 
typical of a surface bubble. In frame 46 the tiny 
central white spot is more diffuse and the enclosing 
semicircular white spots have started to split into 
semicircular lines. In the succeeding frames these 
lines spread out. In frame 54 they are almost con- 
centric circles such as are formed when a pebble is 
dropped in a pool. Some vestige of this configura- 
tion can still be seen in frame 73. 

The fact that the bubble life cycles that have 
been pointed out seem to start after the head of 
water of the drop has vanished must not be con- 
strued as invalidating the second of the explana- 
tions of how cavitation conditions may be realized 
in a liquid drop after impact against a solid surface. 
This was called to my attention by Phillip Eisenberg 
of the Office of Naval Research. It may be neces- 
sary for the negative pressure wave to return to the 
impact surface many times before the bubble nuclei 
have sufficient time to grow to a size that can be 
seen, 
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In evaluating the evidence, which has just been 
presented, for the possible existence of cavitation in 
a waterdrop after collision with a solid surface, it is 
important that the effect of the impact velocity 
should be considered. The terminal velocity of a 
large waterdrop in air is [4] about 26.9 ft/see (820 
cm/sec). One can hardly extrapolate the evidence 
for cavitation found at this velocity to what may be 
found at airplane flight velocities, for example, 
which are of the order of 900 ft/sec. At the high 
impact velocities at which erosion is actually ob- 
served as a result of waterdrop impact the negative 
pressure would be very much increased and_ the 
possibility that cavitation may occur consequently 
much enhanced. This would follow from either of 
the explanations of how cavitation conditions may 
be produced as a result of the impact of a waterdrop 
against a solid surface. 
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FiGgure 9 Schlieren patterns of the collision of a drop of argon saturated water with a glass plate 
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2.4. Time Dependence of the Impact Force 


It is of interest to know something about the ordet 
of time during which the impact force developed by 
the impingement of a freely falling waterdrop 
against a solid surface is effective. A preliminary 
attempt to measure this has recently been made by 
Lawrence Fleming of the Diamond Ordnance Fuze 
Laboratory, using a thin disk of barium titanate 
cera 

The barium titanate disk that was used was about 
lem indiameter. It was coated with silver on both 
sides to make it conducting and was fastened with 
polvstrene adhesive to a small metal base The 
complete unit was then coated with polystyrene 
lacquer to prevent the water of the drop from short 
circuiting the silver lavers. Leads from the silver 
lavers were connected to a Dumont 304 oscilloscope 

The barium titanate unit was placed on a glass 
plate that was located below the exit of the tube 
through which the waterdrops were allowed to fall 
The waterdrops struck the barium titanate disk 
after a free fall of approximately 40 ft. The traces 
that appeared on the oscilloscope screen as a result 
of the waterdrop collisions against the barium 
titanate were photographed directly. One of these 
traces with the time base marker is shown in figure 
10. The time between cycles in the time base marker 


is LOO usec 

It can be seen from the oscilloscope trace that the 
impact force is applied very suddenly. The move- 
ment of the point of light on the escilloscope screen 
due to the increasing force was too rapid to leave a 








Fiaure 10 (xy ‘losco pe trace of the decay of the force produced 
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The decay curve only is 
visible. The force undergoes a rapid decay to zero 


trace in the photograph 
in approximately the space of 1 msec. As has been 
noted in section 2.3, d, this is also approximately the 
time during which the head of water of the drop just 
vanishes, 


2.5. Time Dependence of the Flow Velocity 


The experimental values of radius of flow at 
regular intervals of time for the spread of a drop of 
water after impact at its terminal velocity in-air 
against a glass plate are shown graphically in figure 
11. The data for this graph, which are given in 
table 1, were secured from 50 consecutive pictures 
of the history of the collision of a waterdrop against 
a glass plate. The pictures were the first 50 frames 
of figure 7. The measurements of the diameter of the 
radial flow were made under low magnification with 
dividers and a steel rule graduated to 0.01 in. Each 
experimental diameter is the average of two measure- 
ments. The values plotted on the graph have been 
multiplied by a magnifying factor of 2.7. The 
magnifving factor was determined from a series of 
measurements of the diameter of the maximum 
spread of drops that were obtained from the same 
pipet and that fell approximately the same distance 
The maximum radius of flow found in this way was 
2.3em. The maximum radius of flow in the schlieren 
photographs occurred at frame 44 and was 0.85 ¢m 
The quotient of these values of the maximum spread 
radius produced the magnifying factor. 

Ambiguity exists as to the exact time at which 
the flow was initiated. Any radius of flow that is 
less than the radius of the drop itself is obscured 
by the drop. Extrapolation of the curve indicates 
that flow began not at time equal to zero on the graph 
but at a time equal to approximately 152) usee 
That is, the origin of the curve is translated 152 
usec along the abscissa because, of necessity, the 
radius of flow is zero at zero time, which is taken to 
be the first point of impact. 

The slope of the curve of radius against time is the 
velocity of the radial flow. This curve is also shown 
in figure 1] The observed points are those ob- 
tained directly from the radius-against-time curve 
The velocities are plotted in terms of the terminal 
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velocity, \, of the waterdrop, which was upproxi- 
mately 26.9 ft/see (S820 Expressing the 
radial flow velocity in terms of velocity \° that the 
drop had acquired at the instant at which it impinged 
against the glass plate makes it possible to see at a 
glance how much larger the radial flow velocity is 
than the pact velocity, for times just after the 
collision incident It can also be seen at a glance 
that the radial flow velocity has dropped to the value 
of the impact velocity in less than 1 msec, that is, 
before the head of water that remains of the drop has 
disappeared into the radial flow. To express the 
values of the radial flow velocity in feet per second 
or in centimeters per second it is only necessary to 
multiply the values of the ordinate by 26.9 or by 
S20), respectively. . 

The central and side high lights for this waterdrop 
merged, that is, the head of water of the drop van- 
ished, at the end of 1.2 Furthermore, the 
oscilloscope trace of force as a function of time shown 
in figure 10 for a waterdrop, that was obtained from 
the same pipet and that fell through the same dis- 


Ch sec), 


msec, 
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tance, indicates that the WD pact force vanished at the 
end of | msec. The velocitv-against-time curve 
shown in figure 11 has the same shape as the oscillo 
scope trace of the force as a function of time shown 
in figure 10. This seems to indicate that the radial 
flow is propelled outward by the force of the impact 
\fter the impact vanishes the outward flow 
must continue, although at a much reduced velocity 
under its own momentum. 


foree 


2.6. Collisions of Waterdrops Against Surfaces Hav- 
ing Different Degrees of Resilience or Smoothness 


That the flow characteristics of an 
waterdrop depend strongly on the type of surface 
which the collision occurs as well as on the 


LEER PUTA Se aay 


against 
velocity of the impact can be seen by comparing 
figures 2, 3, and 6 For the study of this effect the 
high-speed camera was placed at a slight angle above 
the surface against which the waterdrop would im 
pinge. The camera was placed rather far from the 
area on the surface where the waterdrop was expected 
to hit in order to increase the probability of getting 
a picture of the collision incident. Unfortunately, 
this resulted in a loss of detail so that enlargements 
of a few frames from the film in the region where the 
collision occurred proved to be of little value The 
films were therefore studied in| motion, using a 
motion-picture projector 


a. Natural Rubber and GR-I Synthetic Rubber 


Impingement of a waterdrop on a sheet of natural 
rubber about SO mils thick and on a sheet of GR-I 
synthetic rubber of the same approximate thickness 
resulted in flow configurations of the same general 
appearance. Quite a bit of spray in upward motion 
was produced in the first stage of the impact. This 
spray moved out with the periphery of the radial 
flow of the waterdrop. See, for example, the com- 
parable configuration that resulted in the case of 
waterdrop impingement against wet filter paper 
shown in figure 3. The upward splash was more 
nearly vertical in the case of waterdrop im pinge- 
ment against natural rubber than in the 
waterdrop unpingement against the S\ nthetic rubber 
Natural rubber is much more resilient than GR-I 
rubber. A steel sphere dropped on a sheet of natural 
rubber undergoes a considerable rebound The same 
steel sphere dropped on a sheet of GR-I synthetic 
| It is possible 


Cuse of 


rubber hardly rises from the surface 
that the difference in resilience of the two rubbers 
may explain the more nearly vertical spray in the 
case of the natural rubber. On the other hand, it 
is possible that the depression formed in the natural 
rubber as a result of the Impact of the waterdrop 
was deeper and that the spray was directed up the 
more nearly vertical sides of the depression to pro- 
duce the configuration of spray that was observed 
The Shore “A” hardness of the natural rubber was 
about 30, whereas that of the GR-I svnthetic rubber 


was about 20 Although this” static hardness 
measurement would seem to indicate that the 
greatest deformation should occur in the GR-I 








taken as a certain 
under dynamic con- 


synthetic rubber, it cannot be 
indication of what to expect 
ditions of loading. 


b. Dry Clean Glass 


The collision and radial flow of a waterdrop on a 
dry clean glass plate was smooth. There was a 
slight spray formation around the head of water 
of the drop as it was driven into radial flow (see fig. 6). 
The glass plates that were used for this study con- 
tained scratch blemishes. 


c. Iron Buffed and Iron as Received 


One side of an iron plate was buffed, the other 
side was retained in the as-received condition. 
Comparison of the flow configurations that resulted 
from the impingement of a waterdrop on each of 
these iron surfaces showed marked differences. In 
the case of the buffed surface there was essentially no 
upwi ard splash of the water as a result of the collision. 
There was only a slight disturbance around the 
vanishing head of water of the drop with a very 
smooth outward radial flow of the water. [mpinge- 
ment of a waterdrop on the iron plate in the as- 
received surface condition resulted in considerable 
upward splash of the water. The amount of spray 
that is formed as a result of waterdrop-to-solid 
collision seems to be a function of the smoothness of 
the surface. The spray formation probably occurs 
at the time of maximum pressure. 

d. Fine Sand in a Petri Dish 


The collision of a waterdrop against fine moist 
sand in a petri dish appeared like an underground 
explosion of dynamite in diminutive. A crown of 
upward moving water splash and sand grains was 
produced by the impact. The print made in the 
sand by the collision consisted of a circular trench 
around a more or less undisturbed mound of sand in 
the center. The central mound of sand probably 
marked the stagnation point of the water flow from 
the drop. The circular trench was dug by the radial 
flow of the water. The configuration is quite inter- 
esting from the standpoint of how impinging water- 
drops may dig pits in soft materials. 


3. Equation for the Impact Pressure Result- 
ing From Waterdrop-to-Solid Collisions 


An early estimate of the mean pressure that re- 
sults from the collision of a solid surface with a water 
sphere was made by Honegger [5]. In his treatment 
the velocity of the center of gravity is the impact 
velocity, V, when the time, f, is zero (at the first point 
of impact). The velocity of the center of gravity 
is |)2 at a time Af later. He assumed that the 
center of gravity during the time Af was displaced 
by the distance d/4, where d is the drop diameter. 
He then applied the impulse momentum equation 
in the form” 


M V/2=Pfat, (1 


where ./ is the mass of the drop, P is the pressure, 


Che notatior ot he ithors has been transcribed into the notat 





and f is the mean contact area between the droy 
and the surface. He obtained in this way the equa- 
tion for the pressure 

P=4X10°°V*, (2 
where P is given in kilograms per square centimeter, 
if V is expressed in centimeters per second. 

De Haller [6] has written as an estimate of the im 
pact pressure when a solid surface strikes a evlinder 
of water from the side the equation for the collision 
of two flat elastic bars. This ts 


P V pres (:3 
, > 
1 + (pi€1, p2C2) 


where p; and p, are the density of the flat water bar 
and of the solid bar, respectively; ¢; and ¢: are the 
speed of sound in water and in the material of the 
solid bar, respectively; and V is the relative speed 
of the water bar with respect to the bar of solid ma- 
terial. This equation neglects the curvature of the 
water cylinder. 

From the results of the experimental work, which 
have been discussed in section 2, it seems reasonable 
to think that a treatment similar to that of the elastic 
impact of solids may provide a reliable estimate of 
the impact pressure. However, the spherical shape 
of the waterdrop, which imposes the condition that 
compressional waves initiated in it by the collision 
are not all started simultaneously, must be taken into 
account, 

Let the following assumptions be made: 

(A) The first effect of the collision ts the initiation 
of a compressional wave that starts to spread through 
the water sphere (see fig. 12). The water molecules 
of the spherical drop in the region traversed by the 
compressional wave attain a velocity in the direction 
in which the impinging solid surface is moving. This 
region of the drop consequently becomes. slightly 
flattened. The velocity taken on by the water mole- 
cules in this region of the drop is dependent on the 
velocity of the impinging solid surface, but it is less 
than this velocity. Only a thin laver of water in 
direct contact with the impinging solid surface takes 
on a velocity that is identical with that of the solid 


Fiat RI 12 Schematic drawing indicating that the water s phe 


is flattened asa result of com pre ssional wavelets initiated in it 
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There is, furthermore, a velocity distribu- 


suriace, 
tion among the water molecules in this region of the 


waterdrop. The following treatment is restricted 
to a consideration of the average velocity given to 
the water molecules in this region. Only the velocity 
in the direction in which the impinging solid surface 
is moving is considered. 

3) Maximum pressure, which is assumed to be 
reached at the end of the time interval Af, is defined 
by the condition that the radiating compressional 
wavelet initiated at a point in the first instant of 
impact should just reinforce the compressional wave- 
let that is started at time Af later in producing mo- 
mentum in the water sphere (see fig. 13). That is, 
the rate of change of momentum of the water mole- 
cules in the region where this reinforcement occurs 
is greater than that of any other of the water mole- 
cules of the drop because it is produced by the com- 
bined effect of the wavelets 

(‘) The second effect of the collision is that the 
first liquid of the drop that encounters the solid sur- 
face during the time interval Af is displaced to form 
a thin evlinder under the drop (see fig. 14). The 
water of the shaded cap of the flattened drop of 
figure 14 is displaced into the ring the triangular cross 
section of which is bounded by the solid surface, the 
vertical boundary of the evlinder of water, and the 
curved boundary of the flattened waterdrop. The 
thin evlinder of water that forms and the region of 
the drop that is traversed by the compressional wave 
during the time interval Af are stippled in figure 14. 
The thin evlinder of water is assumed to be in radial 
flow. The upper surface of it is termed the \-plane. 
The head of water that remains of the drop (and 
which is separated from the surface of the impinging 
solid by the radially flowing thin cvlinder of water) 
retains its shape except for a slight flattening (see 
figs. 2 and 3). 

D) As the collision continues, compressional 
wavelets are initiated at points of contact of the 
remaining water sphere with the A-plane. In this 
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sense the A-plane, or upper boundary of the radial 
flow, may be considered as an effective striking 
surface. 

Consider that the flat surface of a massive solid, 
moving with velocity Vin the 2-direction, Is just at 
the point of collision with a stationary water sphere 
If the first point of contact between the water sphere 
and the flat solid surface is made the origin of a 
rectangular coordinate system, the equation of the 
sphere at the origin of the coordinates is 


+y?+(z2—r’)?=r" (4) 


where the z- and y-coordinates are in the plane of 
the colliding surface, and r’ is the radius of the water 
sphere before the collision. For all points for which 
y=0, eq (4) reduces to the circle 


*\9 ,° ~ 


r+ (2—r')*=r" (5) 


in the r2-plane. The first point of collision between 
the solid surface and the water sphere is given by 
the equation for the circle of contact during the 
collision 


2? +y?—=[R’(2)], (6) 
with /2’(2) equal to zero. 

Consider that the first result of the collision is the 
initiation of a compressional wave that immediately 
begins to spread through the water sphere (see fig. 12 
and assumption A). The water through which the 
compressional wave passes is given a velocity in the 
--direction in which the surface of the massive solid 
was moving. As a result of this, the region of the 
water sphere that has been traversed by the com- 
pressional wave is slightly flattened in this direction. 
The radius of this flattened portion of the water 
sphere is 7, and r is greater than 7’. 

The average velocity of the water in the region of 
the water sphere that has been traversed by the 
compressional wave may be written as aV’, where a@ 


4? 


— 








FIGURI 14. Cross section of the water s phe re that shows sche- 
matically the locus of the vertical components of the compres- 
sional wave and the formation of the thin cylinder of water 


that is in radial flou 








s a coefficient that tells what fraction of velocity V 
is imparted to the water molecules on the average. 
As there is no marked attenuation of the amplitude 
of a compressional wave in water, the coefficient a 
must be governed mainly by the extent of divergence 
of the compressional wave as it spreads through the 
spherical drop As the impact velocity 
for a series of waterdrop-to-solid collisions, the extent 
of flattening in the region of the drop that has been 
traversed by the compressional wave increases, the 
amount of divergence of the wave decreases, and the 
value of the coefficient @ approaches unity 

Consider that the second result of the collision in 
the time interval Af, during which the pressure 
reaches a maximum, is the displacement of the first 
water of the sphere that encountered the solid sur- 
face into a very thin evlinder of water that is in 
radial flow (see fig. 14 and assumption C The 
average velocity at which this displacement occurs 
is(l—a) Vv. It ean readily be shown that in the limit 
the volume of a very thin slice cut off a sphere is 
one-half the volume of a evlinder that has the same 
height and base area. Consequently, in the time 

Af in which the solid surface moving at the displace 
ment velocity (l—a@)Vo has moved a distance 6 
through the water sphere, where 6 is the thickness of 
the cap of the flattened portion of the water sphere 
that is shaded in figure 14, the A-plane has moved a 
4 If the displacement velocity is 


Mncreases 


distance 26 con- 


stant, then the velocity 7 of the \-plane is given by 


Tr 2(1 ~a@ V. 7 


The Massive striking surface IS not al all slowed 
down as a result of the collision It continues to 
move at the velocity Vin the fixed coordinate svs- 


tem In the fixed coordinate svstem, therefore, if 
has traveled the distance 6’. where 
6’= VA R 


through space in the time interval Af. It has, how- 
ever, only displaced the water of the drop which has 


been traversed by the compressional wave and 
which is moving in the same direction that it is 
moving) the distance 

6=(1l—a)VA i) 
Hence 

6 | — a)6’ 1d) 


It can, furthermore, be shown that the radius, +, of 
the flattened portion of the water sphere, through 
which the compressional wave has passed and in 
which the water has the average velocity aV in the 
-direction, is given by 


/ / l a}. 1] 


After the short time Af the solid surface has moved 
a distance 6 through that part of the water sphere 
that has been traversed by the compressional wave, 
the A-plane has moved a distance 26, and the com- 
pressional wavelet that was started at the impact 


instant at point A of figure 13 has moved ¢/F times a 
far as the A-plane has moved, that is, 26¢/7, and has 
reached point B. Here ¢ is the speed of the com- 
pressional wave in water. The radius of the circle 
of contact that is given by eq (6) with ?’(2) replaced 
by R(z) is now sr. This, on assumption B, is the 
point of maximum pressure because at this instant 
water molecules are being accelerated both by their 
own impact in the collision direction 2 and by the 
compressional wavelet that was started at the im- 
pact instant and that has just arrived in_ their 
locality. 

It is possible to solve for 6, the thickness both of 
the shaded cap of water in figure 14 and of the disk 
of water in radial flow in the impact plane after the 
time Af, in terms of the radius, 7, of the flattened 
portion of the water sphere, the velocity of the 
striking surface, V, and the speed of the compres- 


sional wave, c. That is, using the right triangle 
shown in figure 13, 
t¢-6- ™ 

/ —_ o— 15° 1? 

and from eq (5) with 7’ replaced by +, 
r - : r)- 2rz at 13 

so that 

rs $75 — 46 14 


By subtracting eq (12 is seen that 
; | 


6 / . | 5 
c 


According to assumption B, maximum pressure 
after time Af. It is noteworthy that this 
maximum of pressure has been defined by an arbi- 
trary assumption and is not to be regarded aus heces- 
sarily coincident with the maximum of total force in 
figure 10. Actually, since the time Af as evaluated 
in section 4 for a collision of a waterdrop at. its 
terminal velocity in air is 0.04 psec, and since the 
time between cycles of the time base marker in 
figure 10 is 100 ywsee, it is conceivable that it may 
exist at the left of this maximum. Furthermore 
the maximum pressure as it is evaluated at the time 
Af in the following treatment is not an instantaneous 
maximum pressure. To obtain it the total mass of 
water set in motion during the time interval Af, and 
consequently all the force that has acted to produce 
water momentum over the time interval Af, is con- 
sidered The instantaneous force that is acting at 
the last instant of this time interval is not considered 

The average rate of change of momentum of water 
in that portion of the drop that has been traversed 
by the compressional wave over the time interval 
Af is maV Af, where m is the total mass of water in 
the drop that is set into motion as a result of the 
collision. If all of the foree that has acted to 
produce water momentum over the time interval 
Af, at the end of which the radius of the circle of 
contact between the drop and the surface is sy5, 18 
written as mr.,° P?, then 


from eq (14), it 


OcCcUPrs 


[P? mavV lary Al 16 
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is necessary to evaluate the mass of water, m, 
| was set in motion in the 2-direction as a result 
of the collision. After the time Af, part of this 
total mass of water, the mass m’, is moving in hori- 
zoutal flow in the impact plane as a result of the 
impact pressure. The mass moving horizontally is 
given by , 


m ro p(26 I” 957 pd 17 


because in the limit, as has already been noted, a 
slice of thickness 26 and radius «5 cut off a sphere 
is just one-half the volume of a cylinder of the same 
height and base radius. Although the mass of water, 
m’, is moving horizontally after the time Af, it was 
given a velocity in the collision direction 2 during 
this time interval because it was traversed by the 
compressional wave , 
There ts also the mass of water m’’ in the head of 
water remaining of the drop that has a velocity in 
the c-direction because it has been traversed by the 
compressional wave. After the time At the d-plane 
has moved distance 26, and the compressional wave 
has moved to the curved boundary of the stippled 
area Of figure 14. It may be assumed that this 
volume of water has received the signal that collision 
has occurred. The lower boundary of this volume 
of water is the plane 26, where the flattened 
sphere is at the origin of the coordinates. From 
eq (5) the 2-coordinate of the points on the circle of 
radius 7 from which wavelets of compression origin- 
ated ts 
/ / I ; IS 


Az is the distance the A-plane moved up since the 
time that a compressional wavelet was started at the 


corresponding value of 4 Then 
A 26 } } J 19 
Let 2’ be the distance traveled above the plane 26 


Then 


l scav ‘|, [s-r]as. eo 


The locus of the pots VIVES the envelope of the 
components of the compressional wavelets in’ the 
collision direction Substituting the value of A 
Into eq () 


by the compressional wavelet 


from eq (19 


Letting 


and squaring both sides of the equation 


Tr. 442 Bo’ - ; 
z + . — 46° -+- 46) 24 
B D B 


is an expression for 2? in terms of 
volume of water, on substituting this 


Kquation (24 
the enclosed 
value of x*, is 


; ‘ , Ss 7 e - 
7 red Or weet ak O77 |- 20 
a 2 “ ~ — 


The limit of the ratio of this volume over the volume 
of a evlinder of the same height and base area is also 
found to be one-half. Consequently, the mass of 
water m’’ having the average velocity aV in the 
collision direction 2 after the time Af is 


on prp )s 
Mt . d'95 \ 0) 


pw Bars 6. 3(}) 


and the total mass of water that was set into motion 
in the z-direction through having been traversed by 
the compressional wave is from eq (17) and (26 


m—m’+-m"’ =(14-B) 43, 3 pé. (27) 
The pressure, 7’, is then from eq (16) 
P (1-4 B) a3 Tt po| a\ JR) 
I, WAL 
Sitice 
26—2(1—a)VAt=—rAt (29 
and since 
c—v 
re} = 
2 
so that (14+ 8) =e 7, 
» GQ, yy, aie 
/ > i¢ p\ |. ol) 


This is the well-known water-hammer 
multiplied by the factor a@/2, which comes from the 
spherical shape of the waterdrop. It is noteworthy 
that the result that the pressure is some multiple of 
the water hammer pressure because of the spherical 
shape of the waterdrop could have been written down 
at once There would, however, be no rational basis 
for guessing at the value of the multiple. The value 
of the coefficient a found empirically for 
waterdrop-to-solid collisions that can be stopped by 


equation 


ean be 


high-speed motion pictures as is shown in section 4 
For collisions at very high impact velocities, the 
coefficient @ is probably close to units C'onse- 
quently, for these collisions a reliable estimate of the 
impact pressure can also be made from eq (30 
Kquation (30) has verified indirectly by 
de Haller [6], who has measured the impact pressure 
between a flat surface and a jet of water struck from 


been 


the side, using a piezoelectric pressure gage having 


a piston diameter of 1.5 mm. With his expert 
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mental arrangement he was really measuring the 
pressure under short water cylinders struck from the 
side Since the side of a evlinder contains both the 
element of circle of a sphere and the eiement of 
straight line of a plane, the pressure under such 
water shapes should be greater than that given by 
eq (30) but less than that given by eq (3), which is 
the equation de Haller was using for the impact 
pressure. With his piezoelectric pressure gage 


de Haller found the highest pressure value for a 
velocity of 35 m/see to be 310 kg/em?. For this 
velocity eq (3) would predict a pressure of 490 


ke em? if one assumes that the denominator on the 
right-hand side is essentially unity, that the speed 
of sound in water is 5,000 ft/sec, and that V is aV, 
where @ is about 0.9. Equation (30) would predict 
a pressure of 245 kg/cm’, assuming again that a@ is 
about 0.9. The pressure that is developed when a 
flat solid surface strikes a evlindrical water surface 
should be between that which is developed when it 
strikes a flat water surface and that which results 
when it strikes a spherical water surface. De Haller’s 
experimental measurement shows that this is essen- 
tially the case 


4. Equation for the Flow Velocity at! Maxi- 
mum Pressure 


Experimental measurements of the radius of flow 
of a waterdrop at regular intervals of time after 
collision at its terminal velocity in air with a glass 
surface were discussed in section 2.5. Figure 11 con- 
tains the curve of radial flow velocity against time 
The data from which this curve was plotted are given 
in table 1. With eq (30) of section 3 the radial flow 
velocity at the time of maximum pressure can be 
calculated 

Assume that the radial from an impinging 
waterdrop has axial svVmmetry and that the effect of 
Viscosity may Because the pressure 
at the periphery of the radial flow is atmospheric 


flow 
be neglected 


and a be taken aus zero, and because the center of 


the flowing disk is a stagnation point, the equation 


of flow in evlindrical coordinates is 
, pt “" Ov , 
/ vidi | dy M31 
2Jo ‘ Of 
where P is the pressure at the center of the disk. 


of the liquid, / is the observed 


is the obser ed 


p is the density 


velocity of spread, and 7 radius of 


spread. It may be assumed that at the instant of 
maximum pressure Ov/Of is zero. Integration and 
substitution of the expression for P? given by eq 


30) then vields 


acV} @ 5 4 


The value of @ can be calculated from the velocity 
at which the glass surface moved through the water- 
drop. This velocity is (l—a@)V. The central and 
side high lights on the waterdrop for the collision 
incident shown in figure 7 merged in frame 14, 
which was 0.0011, see from the point of impact or 
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time zero. As the average waterdrop from the pipet 
that was used is 0.57 cm in diameter, the velocity at 
which the glass moved through the waterdrop was 
484 cm/sec. 


drop [4] is about 26.9 ft/see (820 em/see Hence, 
l a) S820 184, 
a=0.4. 


From eq (32), assuming that ¢ is about 5,000 ft/sec, 
Ve about 26.9 ft/sec, and using the value of 0.4 for 
a, the radial flow velocity for this waterdrop at the 
time of maximum pressure was 232 ft/sec, or 8.6V. 
In order to plot this point on the velocity-against- 
time curve shown in figure 11 it is necessary to know 
the time at which maximum pressure existed, that 
is, Af must also be evaluated. From eq (9 


6=(1—a) V'At. 
From eq (7) and (15 
6=4r(1—a)"{Vie] 
From eq (7), (9), and (15 
At=4r(1—a) Vie 5 


The radius 7’ of the waterdrop was about 0.29 em, 
so that the radius of the flattened waterdrop from 
eq (11) was about 0.5 em. Using the value of 0.4 
for a, At=0.04 The caleulated velocity is 
plotted on the velocity-against-time curve shown in 
figure 1] and is enclosed in a circle to distinguish it 


uMSeC, 


from the measured values. 

The radial flow of the waterdrop is a direct con- 
of the impact This 
zero in about | msec as is apparent from figure 16 
The head of water of the drop also roes to zero in 
approximately this same time interval, and the 
radial flow velocity drops from 8.6 times the value 
of the impact velocity, V7, to less than WV. After this, 
the water in radial flow continues to move outward 
under its own momentum until the flow is completely 
checked by its internal resistance and by surface 
tension. 


sequence force force “oes to 


5. Equation for the Time Dependence of the 
Flow Radius 


The observed values of the radius of flow at reculat 
intervals of time for the spread of a drop of water 
after collision at its terminal velocity in air against a 
glass plate are shown graphically in figure 11. It is 
possible to caleulate the radius of flow 

At the instant at which it just collides with the 
glass plate, the energy of the falling waterdrop con- 
sists of its kinetic energy and its potential surface 
energy. As soon as if begins to spread on the class 
the kinetic energy begins to be transformed into po- 
tential energy of surface and into dissipated energy 
As long as any part of the head of water remains, the 


The terminal velocity for a large water- | 
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ms of the energy equation are (1) the energy 

netic and potential) remaining in the head of water, 

the potential energy of surface of the radial flow, 

the kinetic energy of the radially flowing water, 
and (4 
time the sum of these energy terms must be equal 
to the total energy of the falling drop at the instant 
of its collision with the glass plate 

rhe first of these energy terms would be difficult 
to evaluate. However, for all times after the head 
of water has vanished, that is, for all times during 
which the radial flow is moving outward under its 
own momentum, this energy term is zero. This is 
the case that will be considered. 

At any instant at which an observation is made, 
the potential energy of surface of the radial flow, 
PE, is given by 


PE=(4 \¢) eI, (34) 


where y, is the surface tension of water against air, y, 
is the interfacial tension between water and glass, 
and ry is the observed value of the radius of spread 

To develop an expression for the kinetic energy it is 
necessary to make some assumption about the shape 
of the cross section of the radial flow. From the 
work of Worthington [1] it appears that the flow of 
a drop of milk is shallowest in the center and deepest 
at the periphery. Figure 15 (a) is an enlarged repro- 
duction of one of the sketches given by Worthington 
As a simple approximation to this shape the assump- 
tion has been made that the radial flow is an extreme- 
ly flat evlinder except for an empty cone that ex- 
tends through half the thickness of the evlinder. 
s shown in 


The assumed cross section of the flow 
figure 15 (b 
The kinetic energy, A“, at the instant when the 


radius is 7) is given by 


where A and ¢ are the thickness of the flow and the 
radial flow velocity, respectively and pis the densit\ 
of the water As vis zero at the center of the flow, 


and ris », at the periphery of the flow, 


W here 


Thickness goes from /,/2 at the center of the flow to 


h, at the periphery, so that 


\/, of the drop, 


In terms of the ass, 


h H.A\/ [Sarprel, 3Y 


the dissipated energy. At any instant of 














(b) 
Figur Ld C'ross section of the flow of al qu dad op 
1) Cro ection of the flow ofa drop of milk as observed by Worthingto b 
sssurmed cro ection of the flow ofa drop of water 


5 (archy Substituting 
and Ay in the expression for 


since the volume of water ts 
the expressions for ¢, h, 
the kinetic energy, 


54 r Mr 1() 
200 °—s«&t* 1 f 

At any instant of observation the dissipated energy, 
Die, is the accumulated dissipation during the time 
interval ¢. This dissipation is 


Dk 2x | Pl Qilt (41) 


where ® is the rate of dissipation of energy per unit 
time per unit volume, and @ is the volume of the 
boundary laver attached to the solid The thickness 
of this laver is ¢, and o is roughly less than ho. 
in an axially symmetrical flow it suffices to consider 
only the variation of the radial flow velocity ¢ with 
the flow thickness 2, for this case . 


since 


Pp Mi ( a ) ’ { $2) 


where w is the VISCOSILA Resorting to dimensional 
analysis and eq (41) and (42), it is seen that 


) reodt, 4°34) 


is the observed radius of flow and C' is a 
constant containing the VISCOSILV pu and having the 
dimensions of viscosity It may be supposed that 
o varies linearly with f, and, since 7, is ro/t, 


where / 


DE k | dt, (44 


where & is a new constant 
The energy equation after the head of water re- 
maining of the drop has vanished is 


MV? a gent 


*t 4 
TT Vet Vea Mr + j- ~° dt 
$f Jol 
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whe re Is the radius of the original drop, and 4 
is the velocity at which the waterdrop strikes the 
surface. Differentiation with respect to time pro- 
aduces 
/) Vireo di My keys 
Dari “ , 4 + () tt) 
dt’ 2t2? dt 2¢ f 
The solution of this equation Is 
) f kK 
¥ 17 
\/ M 
hy T Dr ¥Y , f 
) ) 
where A is an integration constant. The expression 


for the integration constant is found from the con- 
dition that when 7 Pax. t= tmax. Where rg. Indicates 
the maximum spread of the drop and f,,, the cor- 
responding time Substituting the expression for 
the integration constant in eq (47), the equation for 


r, in terms of the time and constants is 


M } 
z + r(y _ > 


The equation is valid for all times after the head 
of water of the drop has vanished until the maximum 
spread is reached It necessary to evaluate the 
constant & from the data For the collision incident 
shown in figure 7 the central and side high lights on 
the waterdrop merge, that the head of water re- 
maining of the drop has vanished at frame 14. At 
frame 14, 1.6, em, f=0.0011, The mass \/ 
of an average drop from the pipet that was used is 


is 


Is, 


Sec 


/ 


0.0989 g¢. The maximum spread was reached at 
frame 44, where » Pear = 2.39 cm and t=t 

0.0040; sec. At present there is no way of measur- 
ing directly the interfacial tension between water 


and glass experimentally, and the value of 22(4 Y 
can only estimated roughly Measurements of 
the surface tension of liquid glass against air have 


he 


been made \lorey 7 lists the following values for 
a commercial soda-lime-silica glass: 404.3 dem at 
1.225° C, 406.2 d/em at 1,125° C, and 407.8 dem at 
1,025° © From these data it appears that the 


value may be about 420 d/em at room temperature 
It is impossible to know and henee to take into ae- 
count any that mav be associated with the 
increased rigidity of the glass on cooling Antonow’s 
rule states [8S] that for two liquids that are mutually 
saturated with each other, the interfacial tension is 
equal to the difference between the surface tensions 
of the two liquids separately To this of 
~2 638d em. This method 


change 


degree 


approximation, 2r(4 


2 
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of approximating the value of (y,.+ y,) Was suggest: 
by William W. Walton of NBS Surface Chemist: 


Section. 


To evaluate *&, eq (48) was put in the form 
M _e were 
272... — 
M M “ 
—k + 2r(4 Y 
2r5 yy a. 


Substitution of the values of AZ, ,,.. ro at frame 14 


fmax. ¢ at frame 14, and 27(4 ye.) produced the 
result. £—0.0076- egs units. 

The values of ro caleulated by use of eq (48) are 
listed in table 1 with the measured values. The 


measured values have been multiplied by a magni 
fying factor. The method of determining this fae- 
tor in 2.5. The calculated 
values are in almost perfect agreement with the em- 
pirical values, which appears to justify, a posteriori, 
the assumptions that were made. The calculated 
are indicated with crosses in figure 11 


discussed section 


is 


values of 7 
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Effect of Darling-Dennison and Fermi Resonance 
on Thermodynamic Functions 


Harold W. Woolley 


Che effect on thermodynamic functions due to resonance of the Darling-Dennison type 
for water and of the Fermi type for carbon dioxide is calculated working directly with 
the elements of the secular determinants that define the energy levels 


1. Introduction 


In calculating thermodynamic functions for various molecular systems, the complication 
is sometimes encountered that the energy levels are expressed indirectly by means of a deter- 
minantal equation. In many cases, the determinant will have nonzero elements only along 
and quite near to the principal diagonal, Frequently, no explicit general solution for the posi- 
tions of the levels as affected by the off-diagonal elements may be known, so that an explicit 
determination of the related thermodynamic properties as influenced by the effect in question 
would seem difficult at best. In some of these cases, useful results can be obtained in a general 
form from the secular equation without making use of explicit individual solutions 


2. Secular Determinant Method 


If the secular determinant, /,,—6,,W), is expanded in polynomial form, beginning with 
the highest power of the unknown, W, then the successive coefficients for the different powers 
of W give the sum of roots and sums of products of the roots, such as 2W,, S2W,W,, SUZW WW, 


i i<j 6< Jl 


ete., when allowance is made for alteration of signs in the successive terms. 
In the calculation of thermodynamic functions, the intermediate step of obtaining parti- 
tion functions involves the evaluation and addition of a large number of Boltzmann factors, 


using the energies of the various levels. This factor, e~” ,/*", is expressible as a series so that 
the summation over a set of energy levels 
(—kT)-™ ' 
ST Op WATS } Ss) oye (1) 
heed ——/ ' hod ' 
i m=) Mt. i 


may be found if the individual sums of powers of energy ZW, 2W?, 2 W%, ete., are known. It 
is possible to find these quantities in terms of sums of products of the roots and hence in 
terms of the coefficients of the polynomial expansion of the secular determinant. In this process 
several algebraic steps are involved including the application of the multinomial theorem 
The result can thus be put in the form 


> W=(> W)?-2 D5 WW (2) 

2S Wis (Qa Wo 3 (By WC WW +3 DED WWW, (3) 
SWCD WAS WH'CHD WywWy 44 cD U V)(U > WWW) 

oe OT ae ee (4) 


SWS W953 WIE WW +5 WHOS WWW 


i ! 1 


5S WDE Wy -5 SW VHHD WWW) 


m 


(DU WW)VVVwWHD) +5 OOH VD WW WWW, (5) 
t k<clom 
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These coefficients for the various terms of the polynomial expansion, > JS SW,W 


SW WwW W,. ete., can be obtained in the Case of ui specific secular determinant and 
i k 


introduced into eq (2) to (6). If the secular equation 
E,,—6,,W|=0 (7 


involves a determinant of the nth order with nonzero elements along the principal diagonal and 
along diagonals immediately adjacent to it, then it is possible to represent the coefficients 
6) so as to represent the quantities YW? directly 


explicitly and introduce them into eq (2) to | 
in terms of the elements of the determinant. Thus for the symmetric determinant with £, 
Rs, E, i+} Fist) ai, kh, {+2 ) w=by By og Feng «=ce,, and with all other elements equal t 
zero, one may with considerable algebra transform from the previous equations to the following: 
>> i > ki (sS 
n—1 
Su > 2S a »S*b Ss 9 
i=] i=l l 
>, Wi=) (64 35a? (ke thy ASS Ok +h +35 07(k, +k 
t i=] i T=1 i=] i=] 
+ 6 4, a 5 ; >> a,b C4 6S a, »b é:. (10 


= ; 
! i 


including terms involving a’s, 6’s and c’s. If the e’s are zero, then 


nm n—l n—2 
= Ws > ki 45 \a?(k? } [- I , f- + 4) b?2(k? } I I j }-? a 
1 i=] i=l i=] 
+ 8S a,a,.,50,(k, +h +h. ..)+2D>5ai- 25° bs 4S ara? ,,+4>5a7b?,, 
l i=! l i=] i 
n—2 n—2 n—s n r ‘ 
45a? b? t 45a, 6? + 15a? ,,b? , >> aa, wie bh. 7 tS) 57d? me (11) 
l i=] l i ti 1 


and if the 6’s also are zero 


So W= D+ 5 a+ Mk, tke tk.) 


and 


, i a I ' vo n—l ; 
DSW D+ 6 Dat dels ot MEE at at Rt BAN BRE + 4k he os + 3K 2) 


' 6 2 lara; i(k? +2k;k ptkk, +2 + 3k? i+2k, + ik; +2t k? 2) 
n—l n—2 n—2 n—3 
‘ \ 6 . X92 . ) ‘2 » 9 > P 
- 25 a°+6> Ja?at,,+6> ata?,,+6>5a?a?,,a?,.. (13) 
l 


Equations (8) to (13) may also be obtained by a matrix method with less effort. This is 
based on a theorem of Borchardt [1]' and of Sylvester [2] which states that if A is a matrix then 
the characteristic roots of A” are the nth powers of the characteristic roots of A. Then if H is 
the energy matrix, the sum of nth powers of roots !V, of the secular equation is given by the 
trace of the nth power of the matrix /7/: 


>> W: =trace H. (14) 


so that eq (1) mas also be written as 


> Wi/kT—trace e~MikT (15) 
and, in effeet, include all the results of eq (8) to (15). 


3. Application of the Method 


Among cases in which molecular energy is given by a secular equation are resonances of the 
Darling-Dennison type and the Fermi type. In the case found by Darling and Dennison [3] 
for the H,O molecule, occurring also in some other cases, the energy W is given by the secular 


equation 


Wwe-W W, 0 0) 

if ws—wW WwW 0 

0) MW Wws—-W Ws 

0) 0 W, Wwe-—Ww ‘ * 0), (16) 


where the energies for the individual interacting levels would have been W; if the perturbation 
had been negligible. The magnitude of the perturbation is indicated by the off-diagonal terms 


Ws iti Wisn, i= Wi 3'ee og42= 2YVE1(1 — 1) 3+ 1) (O34 2). (17) 
Any state indicated by 7, 2, v3 with »,>1 is perturbed appreciably by a state v,—2, v2, v3+2, 
because of the approximate equality of y,; and v3. For higher quantum numbers, large groups 


of levels will be perturbing each other, and the secular equation will involve a secular deter- 
minant of correspondingly high order as indicated by eq (16). In the discussion which follows, 
the number of levels perturbing each other in a single group will be indicated by n, which will 
also be the order of the determinant, The mean value of the energy for the group of levels 
expressed as wave numbers will be indicated as G. The increment in the vibrational partition 
function due to the presence of the perturbation for a group of interacting levels is given by 


' he , l he\? re l he rs 
, kT = . l te 2 = * 3 R 
AQ | pp om Wit5 (Gp) 42 Hi o (gr) 42 we... } (18) 
Fig “ n bracket licate the literature references at the end of this paper 
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According to the formulas derived and the definition of G@, the various increments in the sums 
of powers of the energy W measured from the mean value G will be given by 


A> > W,;=0 
A W!=3 002, (We W paul 


In the remainder of the treatment for this case of resonance, only the effect of A>; W, on 


i 
the partition function and derived quantities will be considered. There is a basis for regarding 
higher order effects, as of A>) W%, ete., as negligible. The levels which are involved in a given 


perturbing set have the same value of »,+ 7; and are subject to the additional restriction that 
either v, or v3 is only odd or only even within a single group. There are thus four cases to con- 
sider: 
I. v, ranges from 2n—2 to 0 while rv; ranges from 0 to 2n—2 
II. », ranges from 2n—1 to 1 while v; ranges from 0 to 2n—2. 
Ill. +, ranges from 2n—2 to 0 while rv; ranges from 1 to 2n—1 
IV. », ranges from 2n—1 to 1 while rv; ranges from 1 to 2n—1 
These four cases give, respectively, for W?, 
I. W: r 2n — 27%) (2n — 274 —1)(27—1) (2%) (20) 
ll. W ; 2n + 1— 27) (2n—27) (2 2) (21) 
lil. W r 2n—27)(2n—2i 27)(2++1) 22 
IV. W ; 2n + 1 —27)(2n— 27) (21) (21+ 1 (23) 
where in each case 7 ranges from 1 ton—1. The evaluationS* W is quite straightforward 
and gives for these cases, respectively, 
y° 2n(n 1) (n | 
*y re tn?—10n+9 24 
) 
») 
. ahivii | fi | _ 
7. T (2n 2n+1 25 
) 
») 
vy” 2n(n—1)(n+1 
_ T: 2n—1)(2n 26 
) 
and 
) 
y° 2nQ0 1) l i. 
y 13 $7) LOn y 2% 
) 


Considering the molecular vibrational energy as if it were due to harmonic oscillators 
having frequencies w, and w; wave numbers, for the two vibrations it is found that G, the average 


energy for the group of n levels is given in the various cases by 


II. G=nw,+(n l)w (29 


Ill. G=(n—1)a,4+-nw (30) 
IV. G=n(@,+ @) (31) 


The expressions are now at hand to permit the estimate of the relative increase in the 
partition function due to the perturbation. Continuing to treat the various cases separately, 
there are obtained for case I, 


h ' 
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Then, in lowest order of approximation, the total effect of the perturbation is indicated by the 
sum of the separate (AQ/Q), for the four cases, to give as the approximate value of A In 
h he 
2(, é r“' ty ‘ oles | ae 
e k1 


AlnQ= 35 (AQ/Q):= 7’ ( >.) 


he h he , , he F he , 
pA g EP 9 TF gg EO 2 9, iF OO 1. gS )| 36) 


If the approximate equality of w, and w; is now taken into account, then letting u=hew/kT 
represent either hew,/kT or hew:/kT or their average value, it is found that 


vy" u-(1 " tied . " ? aye : ; - 
Aln ( i oa ie = Oe Wet pe 4% + De it 1 pe * + thé " = ¢ 
w | f , 
¥Y uné o 
Od 
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This is the approximate contribution to the Gibbs free-energy) function, F°ORT, due to 


the effect of the perturbation. The corresponding contributions to //°/RT and C)/R are 
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7 
u—s-Aln @ and wv’ Aln Q, which give 
u : 


d du 
one _ ‘ —_ ] 1 ¢ 4 
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A value of y of 77.52 em~! has recently been given by Benedict [4] with higher order constants 
and has been used in the present calculation. For the present case, # may be taken as about 
3,650 em~!. As should be expected, the total effect is found to be small. This is due largely 
to the fact that the lowest levels showing this type of resonance are at about 7,300 em7~! so 
that a relatively high temperature is required before the levels in question are appreciably 
excited. The magnitude of the effect is indicated in table 1 for several temperatures. 

One may inquire as to the approximate behavior of the effect in the limiting condition of 
high temperatures. If one neglects higher order effects in the manner of the foregoing deriva- 
tion, whieh would neglect an appreciable relative effect due to anharmonicity, then when the 


functions are expanded in increasing powers of u, the leading term gives the approximate values 
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The quadratic power of temperature is a result of the fourth power dependence of A> W? upon 
the quantum numbers, which supplies a fourth power of 7’ as a factor in the summation. This 
is cut down to the second power by the inverse second power in the coefficient of this term. 

Another case of resonating levels for which details have been worked out is that of Fermi 
resonance in the case of CO,. While the effect of the perturbation upon the thermodynamic 
functions can be expected to be fairly small in this case also, since the perturbation, as in the 
previous case, involves the displacement of some levels upward and some downward in such a 
way as to nearly compensate, it is nevertheless a larger effect than in the case of H,O, mainly 
because the energy of the lowest resonating levels is near 1330 em~' for CO, as compared with 
the 7,300 em~ for H,O, so that comparable excitation takes place at a lower temperature. 

It was shown by Dennison [5] and Adel and Dennison [6] that a satisfactory representation 
of the CO, spectrum could be obtained if account were taken of perturbations within groups of 
levels with equal values of /, the azimuthal quantum number, and with equal values of 27, + 2. 
Then the effect of the perturbation within a group of n levels of this kind is indicated by the 
shifts of roots of such a secular equation as (16) with corresponding effects on the thermodynamic 
functions indicated in terms of the partition function according to eq (18) and (19). In each 
group of n levels, the value of v7, ranges from zero to n—1, with 7 going from its maximum value 
to its minimum value which is equal to /|. The number of levels interacting within one group 
is given by n=3(V—(/ 1. The elements of the determinant which are adjacent to the main 
diagonal were presented by Dennison in a form which can be so converted as to indicate that 


W? ..,= W2i(n—1)(V+2—n—-12), (42) 


where 7 ranges from 1 to n—1. For the evaluation of the principal contribution to the parti- 
tion function according to eq (18), there is needed a knowledge of the distribution of the 
unperturbed vibrational levels. It is assumed that these are given in the somewhat uncon- 


ventional form 
G,= Soest SX, (e;— 1) + DO X yw). (43) 


Then the average value of G for a given group of levels is given by 
G=w.V+ w3r3+ w, (n— 1) + .X2(V?7— V) + X93(¢3—v3) +Xan(n—1)° 
Xo.V0r.+ Xy,V(n—1)+X,3(n—1)v,, (44) 
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By simple summation, 
SW: w2 a Ee 2)—5n | (46) 
and for the next term of (18) 
STW? (Wet We, 2 ake i) oo (47) 











where 6,=. X,,—5 Ve—. Ny. After multiplying by the exponential Boltzmann factor for the 


3 3 6 

average Vibrational energy, the summation over Vand x» follow, but with the requirement that 
half values of the summand are to be used when V' has its lowest value, 2(m—1), in order to take 
account of the absence of / degeneracy for /=0. The choice of doing this rather than multiply- 
ing by 2 for /#0 comes about as the result of the fact that half of the rotational levels are 
absent due to the symmetry of the molecule and the zero nuclear spin for oxygen. Then this 
is compared with a total approximate Q similarly omitting any factor for / doubling. If the 
effects for anharmonicity are disregarded, the leading term is obtained as 


+S 
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The effect of Fermi resonance, using u=hew/kT and disregarding the distinction between 
and 2w., is thus given approximately [7] by 
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The effect can be obtained in more detail by including an additional contribution from eq 


46 Using u;=hew, kT and using w,—2a, after the first term, eq (48) is replaced by 
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A calculation was made for CO, using as values for the constants (in em™') 


w, = 1342.86 w) = 666.27 w3 = 2349.15 
An 2.20 Ao 0.75 X33 12.46 
NX }2=3.76 Ay 18.54 No, 12.48 
with B+ g.= 1.42 and Wo=51.01 em, giving w,=3.39, X,, 4.24, X,,=3.38, and X,;=3.21. 


With these constants, the magnitude of the effect as given by eq (48) is indicated in table 2 as 
the first column under each heading while the combined effect consistent with eq (52) is shown 
as the second column under each of the headings, A(—F°/RT), A(H°/RT) and AC?/R. The 
entire effect for CO, is rather small and in many computations could be neglected. For elevated 
temperature, the leading term in this case is roughly linear in the absolute temperature rather 
than quadratic as was found to be the case for H,O. This is due to the essentially cubic 
dependence of W?,., upon quantum numbers in the present case, as compared with the 
quartic dependence in the case of H,O. It is also seen that this simple dependence upon tem- 
perature is modified greatly by the effect of anharmonicity. The large effect of anharmonicity 
at high temperature can be regarded as due to its lowering of entire groups of resonating Jevels 
with a resultant increased population in such groups, which provides an increase in thermody- 
namic effect. 

This treatment does not make any adjustment for the effect of resonance on the composi- 
tions of states, since the rotational perturbation is quite small, and, in addition, the resonance 
shifts have compensating effects in first approximation as has been noted in deriving the larger 
contributions due to perturbation which have been taken into account. They would thus be 
quite negligible 


TABLE 2 Effect of Fermi resonance on thermodynamic functions for CO, 
remperatu (FRI HIRI MC yi 
A 
2K 5x10 36x10 27X10 
2 22 10 1210 70X10 
si 5x10 2X10 llixlo-* 
1K Is 10 X10 17x 
On 33x 10 85x10 295 10-4 
(49 (92 1(50 eq (52 1(51 (92 
1M) 0.0012 0. 0012 0. OO1S 0. OOLS 0. OO33 0. OOBS 
2 000) 0029 0030 0033 0036 O82 OUTS 
(MK) O45 OO4S OO4s 0056 OO9S Olly 
in) OO8! O87 0063 OOT7S 0124 O173 
LD (W77 , OST OOTS O103 OL55 0233 
Equations (49), (50), and (51) also represent approximately the effect of Fermi resonance 
in H.O involving levels for which 27;+ 7, is a constant, even though v, is not degenerate as for 
CO,. Benedict [4] gives the matrix element (7, 72+-2)7)+1, 72) = 4! 200(7) + 1)(v2+2)(+1) }4, 
indicating that W2—200 em~*. A few values for the Fermi resonance effect for H,O are given 


intabie 3. The effect is much smaller than the Darling-Dennison resonance effect at the higher 


TARLE 3 Effect of Fermi resonance on thermodynamic functions for HO 
I peratu \ F°/RT \CIT?| RT AC I 
A 

AM) xx 10°" 10 10 

OOo +x 10 kx 10 2x10 

2) 10 2x10 x10 

4,000) 2x If x10 5x10 

1,000 4x10 ™ 10 7x10 

inn) ix 10 ™ 10 ux 10 
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temperatures although it is the larger at 1,000° K and below. Its earlier rise is partly due to 
the resonance occurring first near 3,400 em~! as compared with 7,300 cm~! for the Darling- 
Dennison resonance. The magnitude of the resonance effects for H2,O and CO, are shown in 
figures 1 and 2, respectively. 
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FIGURE ] Effect of resonance on thermodynamic Figure 2 Effect of Fermi resonance on thermody- 
functions for H,O namic functions for CO 
hi, fi: Contribution of Darling-Dennison resonance to C,/R ce, h, f: Contributions to C./R, (iH I RT, and I I R1 
H°—Ff RT, and F°-F RT ce, h, f: Contribution of Upper branches with anharmonicity, lower branches neglectit 
Fermi resonance to ¢ R, (il I RT, and } I RT inharmonicit ys 


4. Conclusion 


A useful method of calculation of thermodynamic effects for levels defined by a secular 
determinant equation has been demonstrated 

Many other instances with similar perturbations are known and their effects on thermo- 
dynamic functions should be calculable by similar methods if the theoretical description of the 
perturbation has been worked out and the spectra fitted in the detail that has been achieved 
for H,O and CO,. In particular, C,H, gives an example of Darling-Dennison resonance for 
which eq (37), (38), and (39) would apply. Other molecules with Fermi resonance similar to 
that for CO, to which eq (49), (50) and (51) apply include Cs,, C,H,, C.D, and N,O 
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